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Abstract 

Sea  ice  is  a  transiucent  materiai  with  an  intricate  structure  and  compiex 
apticai  properties.  Understanding  the  retiection,  absorption,  and  transmis¬ 
sion  ot  shortwave  radiation  by  sea  ice  is  important  to  a  diverse  array  ot 
scientitic  probiems,  inciuding  those  in  ice  thermodynamics  and  poiar  ciima- 
toiogy.  Radiative  transter  in  sea  ice  is  a  combination  ot  absorption  and 
scattering.  Ditferences  in  the  magnitude  ot  sea  ice  opticoi  properties  are  due 
primariiy  to  ditferences  in  scattering.  Spectrai  variations  are  mainiy  a  resuit  of 
absorption.  Changes  in  such  opticai  properties  as  the  aibedo,  refiectance, 
transmittance,  and  extinction  coefficient  are  directiy  reiated  to  changes  in  the 
state  and  structure  of  the  ice.  Physicai  changes  that  enhance  scattering,  such 
as  the  formation  of  air  bubbies  due  to  brine  drainage,  resuit  in  iarger  aibedos 
and  extinction  coefficients.  The  aibedo  is  quite  sensitive  to  the  surface  state,  if 
the  ice  has  a  snow  cover,  aibedos  are  iarge.  in  contrast,  the  presence  of 
iiquid  water  on  a  bare  ice  surface  causes  a  decrease  in  aibedo,  which  is 
more  pronounced  at  ionger  waveiengths.  Sea-ice  opticai  properties  depend 
on  the  voiume  of  brine  and  air  and  on  how  the  brine  and  air  are  distributed. 


Cover:  Dr.  T.C.  Grenfell  measuring  melt  pond  albedos  on  first-year  Ice 
near  Barrow,  Alaska.  A  Kipp  radlametep  which  measures  fatal 
shortwave  Irradlance,  Is  In  the  faregraund.  The  cylindrical  Instru¬ 
ment  an  the  trIpad  Is  a  scanning  spectraradlameter  that  measures 
spectral  Irradlance  tram  400  to  2500  nm. 
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The  Optical  Properties  of  Sea  Ice 


DONALD  K.  PEROVICH 


INTRODUCTION 

Sea  ice  is  a  translucent  material  with  an  intri¬ 
cate  structure  and  complex  optical  properties.  Un¬ 
derstanding  the  reflection,  absorption,  and  trans¬ 
mission  of  shorfwave  radiafion  by  sea  ice  is 
imporfanf  fo  a  diverse  array  of  scienfific  prob¬ 
lems.  If  is  of  fundamenfal  concern  in  freafing  large- 
scale  problems  in  ice  fhermodynamics  and  polar 
climafology.  The  summer  melf  cycle  of  fhe  Arcfic 
sea  ice  cover  is  driven  by  shorfwave  radiafion, 
making  fhe  inferacfion  of  shorfwave  radiafion 
wifh  sea  ice  a  crifical  componenf  of  fhe  heaf  bal¬ 
ance  of  fhe  ice  cover  (Maykuf  and  Unfersfeiner 
1971,  Maykuf  and  Perovich  1987,  Thorndike  1992, 
Eberf  and  Curry  1993).  Of  particular  climatologi¬ 
cal  concern  is  undersfanding  fhe  sea  ice  albedo 
feedback  mechanism  (Ingram  ef  al.  1989).  During 
fhe  summer  fhe  ice  cover  begins  fo  melf  due  fo 
fhe  inpuf  of  solar  radiafion.  This  melfing  fends  fo 
decrease  fhe  surface  albedo  and  increase  fhe  heaf 
inpuf,  fhereby  accelerating  fhe  melf  process.  Be¬ 
cause  of  fhe  climatological  inferesf  in  fhe  heaf 
balance  of  sea  ice,  fhere  is  also  a  need  for  large- 
scale  spatial  and  femporal  information  on  ice  pack 
albedos.  Properly  inferprefed,  fhe  reflected  radi¬ 
ance  measured  by  visible  and  near-infrared  safel- 
life  sensors  can  provide  such  informafion.  In  ad¬ 
dition,  fhe  amounf  and  specfral  composition  of 
shorfwave  radiafion  fransmiffed  fhrough  sea  ice 
sfrongly  impacfs  primary  producfivify  and  bio¬ 
logical  acfivify  in  and  under  a  sea  ice  cover  (Soo 
Hoo  ef  al.  1987,  Arrigo  ef  al.  1993).  Visible  lighf 
benefifs  ice  biofa  by  confribufing  fo  phofos5mfhe- 
sis,  white  ulfraviolef  lighf  can  damage  organisms. 

This  monograph  focuses  on  fhe  opfical  proper¬ 
ties  of  sea  ice.  The  goal  is  fo  provide  an  infroduc- 
fory  fuforial  fo  fhe  topic,  nof  fo  be  a  complefe 
compendium  of  work  in  fhe  field.  The  physical 
principles  underlying  radiative  fransfer  in  sea  ice, 
including  scattering  and  absorption,  are  discussed. 


along  wifh  fhe  imporfance  fo  opfics  of  fhe  sea  ice 
physical  sfafe  and  sfrucfure.  Observational  resulfs 
are  presenfed,  wifh  fhe  emphasis  placed  on  ex¬ 
plaining  fhe  wide  variabilify  in  sea  ice  opfical 
properties  in  ferms  of  ice  physical  properties  and 
radiative  fransfer  fheory.  An  overview  is  given  of 
exisfing  sea  ice  radiafive  fransfer  models  presenf- 
ing  fheir  basic  characferisfics,  solufion  schemes, 
sfrengfhs,  and  limifafions.  Finally,  currenf  research 
areas  and  problems  of  inferesf  in  sea  ice  opfical 
properties  are  discussed.  Since  fhe  presence  of  a 
snow  cover  can  greafly  impacf  lighf  reflecfion  and 
fransmission  fhrough  sea  ice,  some  menfion  is 
made  of  fhe  opfical  properties  of  snow.  An  excel- 
lenf  review  of  fhe  opfical  properties  of  snow  is 
provided  by  Warren  (1982).  The  opfical  proper¬ 
ties  of  ice  biofa  and  parficulafes  found  in  fhe  ice 
(Arrigo  ef  al.  1991,  Roesler  and  Ifurriaga  1994) 
are  also  discussed  briefly  because  of  fheir  impacf 
on  radiafive  fransfer  in  sea  ice. 


BACKGROUND 

By  "opfical"  we  refer  fo  fhe  porfion  of  fhe  elec- 
fromagnefic  specfrum  fhaf  is  coincidenf  wifh  fhe 
wavelengfh  range  of  radiafion  from  fhe  sun,  from 
roughly  250  nm  fo  2500  nm  (Fig.  1).  The  solar 
porfion  of  fhe  elecfromagnefic  specfrum  is  also 
referred  fo  as  shorfwave  radiafion.  The  opfical 
region  can  be  divided  info  fhree  segmenfs:  ulfra- 
violef  lighf  from  250  fo  400  nm,  visible  lighf  from 
400  fo  750  nm,  and  near-infrared  lighf  from  750  fo 
2500  nm.  The  ulfraviolef  can  be  furfher  divided 
info  UV-C  from  200  fo  280  nm,  UV-B  from  280  fo 
320  nm,  and  UV-A  from  320  fo  400  nm.  Because  of 
sfrong  absorption  in  fhe  afmosphere,  essenfially 
no  UV-C  reaches  fhe  Earfh's  surface.  If  is  in  fhe 
UV-B  where  lighf  levels  are  subsfanfially  enhanced 
by  fhe  depletion  of  sfrafospheric  ozone  (Frederick 
and  Lubin  1988,  Lubin  ef  al.  1989,  Tsay  and 


Stamnes  1992,  and  Smith  et  al.  1992a)  and  can 
have  a  deleterious  impact  on  living  organisms 
(Smith  1989,  Smith  et  al.  1992b).  The  familiar  spec¬ 
trum  of  visible  lighf  is  also  shown  in  Figure  1 
from  violef  (400  nm)  fo  blue  (450  nm)  fo  green 
(550  nm)  fo  yellow  (600  nm)  fo  red  (650  nm). 

"Properfies"  refers  fo  fhe  paramefers  fhaf  are 
used  fo  describe  fhe  reflection,  absorption  and 
fransmission  of  solar  radiation  by  sea  ice.  The 
ferminology  of  radiafive  fransfer  is  infricafe  and 
voluminous.  If  also  has  fhe  unforfunafe  affribufe 
fhaf  fhe  same  physical  quantify  may  have  a  dif- 
ferenf  name,  depending  on  whefher  an  oceanog¬ 
rapher,  an  asfrophysicisf  or  a  biologisf  is  speak¬ 
ing.  To  avoid  a  Babel  of  jargon  we  shall  limif 
ourselves  fo  fhe  ferms  needed  for  a  basic  under- 
sfanding  of  fhe  optical  properfies  and  shall  fol¬ 
low  fhe  ferminology  conventions  of  fhe  sea  ice 
liferafure. 

The  specfral  radiance  f(0,(|),A,)  is  fhe  power  in  a 
ray  of  lighf  in  a  parficular  direcfion,  where  0  is 
fhe  zenifh  angle  (0  pointing  downward,  n  poinf- 
ing  upward),  (|)  is  fhe  azimufh  angle  and  X  is  fhe 
wavelengfh.  The  specfral  radiance  is  defined  as 
fhe  radianf  flux/ nanomefer  per  unif  area  per  unif 
solid  angle  in  a  parficular  direcfion  and  has  unifs 
of  W  m“^  sr“^  nm“^.  The  specfral  irradiance  F{X)  is 
simply  fhe  radiance  projecfed  onfo  a  plane  sur¬ 
face  and  infegrafed  over  a  hemisphere.  Because 
of  fhis  projecfion  fhe  radiance  is  scaled  by  cos  0. 
The  downwelling  irradiance  the  radiance 

infegrafed  over  downward  direcfions  (e.g.,  from 
fhe  sky),  and  fhe  upwelling  irradiance  F„(A,)  is  fhe 
radiance  infegrafed  over  upward  direcfions  (e.g., 
from  fhe  surface).  This  can  be  expressed  formally 
as: 


2jt  7i/2 

^d(^)=  I  I  f  (0,(|),A,)cos0  sin0d0d|) 

6=0  0=0 

271  71 

^u(^)=  I  I  f(0,(|),A,)cos0sin0d0d|)  • 

6=0  e=7t/2 

The  mosf  sfudied,  and  mosf  used,  optical  prop¬ 
erly  of  sea  ice  is  fhe  albedo  (a).  The  specfral  al¬ 
bedo  is  simply  defined  as  fhe  fracfion  of  fhe  inci- 
denf  irradiance  fhaf  is  reflecfed: 

^d(0,?L) 

where  fhe  0  designafes  fhe  surface.  In  sea  ice  fher- 
modynamic  sfudies  fhe  wavelengfh-infegrafed,  or 


fofal,  albedo  is  offen  a  quanfify  of  inferesf, 
since  if  is  a  measure  of  fhe  fofal  solar  energy 
absorbed  by  fhe  ice  and  ocean  (Maykuf  and 
Unfersfeiner  1971,  Maykuf  and  Perovich  1987).  If 
can  be  expressed  in  ferms  of  fhe  specfral  albedo 
and  fhe  specfral  incidenf  irradiance  as 


|a(A,)Fd(0,A,)dA, 
I  F^{0,X)dX 


(1) 


The  fofal  albedo  depends  on  fhe  specfral  disfri- 
bufion  of  fhe  incidenf  irradiance  as  well  as  on  fhe 
specfral  albedo  of  fhe  surface.  Thus  a  change  in 
cloud  conditions,  and  fhereby  fhe  incidenf  spec¬ 
fral  irradiance,  can  resulf  in  changes  in  fhe  fofal 
albedo  (Grenfell  and  Maykuf  1977). 

For  some  problems  a  knowledge  of  fhe  angu¬ 
lar  disfribufion  of  fhe  reflecfed  radiance  is  needed. 
For  example,  in  climafe  sfudies  if  would  be  use¬ 
ful  fo  derive  large-scale  ice  albedos  from  safellife 
dafa.  However,  safellife  sensors  have  narrow  fields 
of  view  and  measure  reflecfed  radiance.  The  key 
fhen  is  fo  relafe  fhe  radiance  reflecfed  af  fhe  view¬ 
ing  angle  of  fhe  insfrumenf  fo  fhe  albedo  of  fhe 
ice.  In  order  fo  do  fhis  fhe  angular  disfribufion  of 
reflecfed  radiance,  characferized  by  fhe  bidirec- 
fional  reflecfance  disfribufion  function  (BRDF), 
musf  be  known.  The  formal  definition  of  fhe  BRDF 
is  (Nicodemus  ef  al.  1977,  Warren  1982,  Perovich 
1994) 


dl  (0,(|),A,) 

cos(0o)dF(0o,(|)o,A.) 


where  0o  and  are  fhe  solar  zenifh  and  azimufh 
angle,  F(0q,  (|)q.  A.)  is  fhe  incidenf  specfral  irradi¬ 
ance,  and  R  has  unifs  of  sferadians”^.  Formally  R 
is  a  derivafive  quanfify,  similar  fo  a  probabilify 
densify  funcfion,  defined  in  ferms  of  infinifesi- 
mal  angles.  In  practice,  fhe  definifion  is  exf ended 
fo  finife,  measurable  angles,  so  fhaf  dI—>AI  and 
dF^AF. 

Lighf  fransmission  fhrough  fhe  ice  is  charac¬ 
ferized  by  fhe  fransmiffance  T(A),  which  is  simi¬ 
lar  fo  fhe  albedo  in  fhaf  if  is  fhe  fracfion  of  fhe 
incidenf  irradiance  fhaf  is  fransmiffed  fhrough 
fhe  ice.  Lighf  affenuafion  in  fhe  ice  is  offen  repre- 
senfed  using  an  irradiance  exfincfion  coefficienf 


k(z,A) 


-1  dF^{z,X) 
F^  (z,  a)  dz 


where  F^{z,  A)  is  fhe  downwelling  specfral  irradi¬ 
ance  af  depfh  z  in  fhe  ice. 

Lef  us  now  examine  fhe  difficulfies  in  defer- 


2 


I  pn 


I 


r  -u 


I  itiiaic-UHih  I 


H 


I 

in 


itcnn 


icpu  fw-rwQr 

-Mm 

LVA  SD-UUn 

Lw  B  an  ■ )»« 


Figure  1.  The  optical  portion  of  the  electromagnetic  spectrum.  Visible  light  is 
from  400  nm  (violet)  to  750  nm  (red). 


mining  the  optical  properties  of 
sea  ice  by  posing  a  simply  stated 
question:  "What  is  the  albedo  of 
sea  ice?"  Albedos  are  straight¬ 
forward  to  determine.  A  radi¬ 
ometer  is  used  to  measure  the 
irradiance  incident  on  a  surface 
and  reflected  from  the  surface. 

The  albedo  is  constrained  to  lie 
between  0,  if  none  of  the  inci¬ 
dent  irradiance  is  reflected,  and 
1,  if  all  the  incident  is  reflected. 

At  first  glance  this  appears  to  be 
an  easy  question  to  answer. 

Figure  2  is  an  aerial  photo¬ 
graph  of  a  small,  roughly  one- 
quarter-square-kilometer,  area  of 
a  typical  summer  Arctic  scene.  The  melt  season 
has  begun  and  there  is  a  tremendous  amount  of 
spatial  variability  in  ice  surface  conditions:  snow- 
covered  ice,  bare  white  ice,  blue  melt  ponds, 
dirty  ice,  and  areas  of  open  water.  This  vari¬ 
ability  is  also  manifested  in  the  wavelength- 
integrated  albedo,  which  ranges  from  0.05  for 
open  water,  to  0.2  to  0.4  for  ponded  ice,  to  0.5 
to  0.7  for  bare  ice,  to  0.75  to  0.85  for  snow- 
covered  ice.  Observations  of  wavelength-inte¬ 
grated  albedo  for  a  full  range  of  sea  ice  types 
and  conditions  are  summarized  in  Figure  3. 
Considerable  variability  in  albedo  is  apparent. 
Determining  that  the  albedo  falls  between  0.05 
and  0.9  still  does  not  provide  an  adequate  an¬ 
swer  to  our  question  of  "What  is  the  albedo  of 
sea  ice?"  Indeed,  while  the  question  is  simple 
to  state,  it  is  extremely  difficult  to  answer  on  a 
large  scale. 

Considering  the  complicated  and  variable 
physical  structure  of  sea  ice,  variability  in  the 
optical  properties  should  not  be  surprising.  To 
understand  and  explain  this  variability,  it  is  neces¬ 
sary  to  examine  the  physical  state  and  structure 

Bulk  Albedo 


Figure  2.  Aerial  photograph  of  typical  Arctic  summer 
scene  taken  from  an  altitude  of  600  m  on  3  August 
1994  at  78°N,  177°W.  The  horizontal  extent  is  approx¬ 
imately  425  m. 

of  the  ice.  As  Figure  2  indicates,  sea  ice  exhibits  a 
great  degree  of  horizontal  variability  with  diverse 
surface  conditions,  including  ponds,  bare  ice,  and 
snow-covered  ice,  and  thicknesses  that  range  from 
open  water  to  pressure  ridges  over  10  m  thick. 

There  is  also  vertical  complexity, 
with  ice  properties  such  as  tem¬ 
perature,  salinity,  brine  volume 
and  air  volume  changing  signifi¬ 
cantly  from  the  ice  surface  to  the 
ice/water  interface.  The  details  of 
sea  ice  physical  properties  and 
structure  are  summarized  in 
Weeks  and  Ackley  (1982).  What  is 
most  germane  to  optics  is  that  sea 
ice  has  an  intricate  structure  con¬ 
sisting  of  an  ice  matrix  with  in¬ 
clusions  of  air,  brine,  solid  salts 
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Figure  3.  Range  of  observed  values  of  total  albedo  for  sea  ice.  The  albedos 
are  from  Burt  (1954),  Chernigovskiy  (1963),  Langleben  (1971),  Grenfell 
and  Maykut  (1977),  and  Grenfell  and  Perovich  (1984). 
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and  contaminants,  and  that  it  is  a  material  that 
exists  at  its  salinity-determined  melting  point. 
Therefore,  changes  in  temperature  result  in 
changes  in  its  physical  properties  and  structure. 
One  of  fhe  goals  of  fhis  fuforial  is  fo  illusfrafe,  af 
leasf  qualifafively,  how  changes  in  fhe  ice  physi¬ 
cal  properfies  are  relafed  fo  changes  in  opfical 
properfies.  To  accomplish  fhis,  we  musf  firsf  ex¬ 
amine  fhe  fheorefical  underpinnings  of  radiafive 
fransfer  in  sea  ice. 


THEORY 

The  inferacfion  of  solar  radiafion  wifh  sea  ice 
is  illusfrafed  schemafically  in  Figure  4.  The  inci- 
denf  radiafion  field  consisfs  of  a  direcf  beam  com- 
ponenf  from  fhe  sun  and  a  diffuse  componenf 
from  fhe  sky  and  clouds.  If  if  is  complefely  cloudy 
and  fhe  solar  disk  is  nof  visible,  fhe  incidenf  ra¬ 
diafion  field  is  considered  fo  be  diffuse.  Depend¬ 
ing  on  sky  and  surface  condifions  some  porfion 
of  fhe  incidenf  radiafion  is  specularly  reflecfed 
from  fhe  surface.  A  porfion  of  fhe  incidenf  radia¬ 
fion  is  reflecfed  from  fhe  ice,  a  porfion  absorbed 
in  fhe  ice,  and  a  porfion  fransmiffed  fhrough  fhe 
ice.  As  we  shall  see,  fhe  relafive  sizes  of  fhese 
porfions  are  dependenf  on  fhe  physical  proper¬ 
fies  of  fhe  ice  and  on  fhe  wavelengfh  of  fhe  lighf. 

Af  opfical  wavelengfhs,  radiafive  fransfer  in 
sea  ice  is  governed  by  fwo  processes:  absorpfion 
and  scaffering.  As  a  ray  of  lighf  passes  fhrough 
sea  ice,  some  of  fhe  lighf  is  absorbed  by  fhe  ice 
and  some  of  if  is  scaffered  from  fhe  beam  in  dif- 
ferenf  direcf  ions.  This  is  expressed  more  formally 


Figure  4.  Schematic  of  radiative  transfer  in  sea  ice. 


as  fhe  equafion  of  radiafive  fransfer  for  a  plane 
parallel  medium  (Chandrasekhar  1960): 

-p — ^  =  J(i:,p,(|),A,)-S(i:,p,(|),A,)  (2) 

where  I  =  fhe  radiance 

p  =  fhe  cosine  of  fhe  zenifh  angle  6 
(|)  =  fhe  azimufh  angle. 


Scaffering  is  included  in  fhe  S  ferm,  which  is  re¬ 
ferred  fo  as  fhe  source  function,  x  is  fhe  nondi- 
mensional  opfical  depfh  and  is  defined  as 

T  (A.)  =  [k  (A.) -I- a  (A,)]  z 


where  k  is  fhe  absorpfion  coefficienf,  cr  is  fhe  scaf¬ 
fering  coefficienf,  and  z  is  fhe  physical  depfh.  The 
single  scaffering  albedo 
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A  (A) -I- a  (A) 


gives  fhe  fracfional  loss  due  fo  scaffering  (Chand¬ 
rasekhar  1960,  Mobley  1994).  atg  ranges  from  0  for 
a  purely  absorbing  medium  fo  1  for  a  purely  scaf¬ 
fering  medium.  A  plane  parallel  medium  is  hori- 
zonfally  homogeneous,  buf  can  have  verfical 
variafions. 

The  compacf  form  of  eq  2  belies  ifs  frue  com¬ 
plexify.  This  complexify  becomes  evidenf  if  fhere 
is  scaffering  in  fhe  medium  (ajg  >  0)  and  fhe  source 
function  is  expressed  in  defail.  For  a  plane-paral¬ 
lel  medium  wifh  a  direcf  incidenf  beam,  fhe  source 
function  is  expressed  as 


S(t,p,(|),A) 


G5  ^ 

^  I  I  P  (P.  Ph  ¥)  I  ('c,  P,  A) 

_i  0 


d^'d(^-  p (p, ph (|), (^') e 

where  p(p,p',(|),(|)')  is  fhe  phase  funcfion  and  Eg  is 
fhe  radiance  of  fhe  direcf  beam  componenf  of  fhe 
incidenf  radiafion  field.  Wifh  scaffering  included, 
eq  2  is  an  infegro-differenfial  equafion  and  is  nof 
readily  amenable  fo  solufion.  However,  while  if 
is  difficulf  fo  solve  fhe  equafion,  if  is  sfill  sfraighf- 
forward  fo  undersfand  qualifafively.  The  double 
infegral  ferm  is  used  for  diffuse  radiafive  pro¬ 
cesses  only  and  represenfs  scaffering  of  fhe  radi¬ 
ance  field  J(t,pA(|)',A)  from  differenf  direcfions  info 
fhe  direcfion  of  fhe  solufion  (p,0).  How  much  of 
fhis  lighf  is  scaffered  from  one  direcfion  fo  an- 
ofher  is  defined  by  fhe  phase  funcfion  p(p,pA(|),(|)'). 
The  phase  funcfion  is  normalized  so  fhaf  ifs  in- 
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tegral  over  angle  is  equal  to  one.  The  second  term 
provides  the  contribution  of  scattered  light  from 
fhe  affenuafed  direcf  beam  Eq(X).  This  ferm  is 
needed  only  if  fhere  is  direcf  incidenf  irradiance 
as  well  as  a  diffuse  componenf. 


Absorption 

It  is  time  to  examine  absorption  and  scattering 
in  sea  ice  in  detail,  starting  with  absorption  be¬ 
cause  it  is  the  simpler  of  fhe  fwo  processes.  Con¬ 
sider  fhe  case  of  absorpfion  only  for  a  direcf  beam 
of  lighf  normally  incidenf  on  a  medium.  Since 
fhere  is  no  scaffering,  S  =  0  and  a  =  0.  For  normal 
incidence,  6  =  0,  which  means  p  =  1.  Equafion  2 
fhen  reduces  fo 


-dl  {z,X) 
k{X)dz 


l{z,X) 


which,  when  solved,  gives  fhe  familiar  exponen¬ 
tial  decay  law 

l{z,X)  =  l{0,X)e-^^  (3) 

also  known  as  Beer's  law  or  fhe  Bouguer-Lam- 
berf  law.  Radiative  fransfer  in  a  purely  absorbing 
medium  is  quife  simple  fo  describe.  The  radiance 
decreases  exponenfially  wifh  depfh  in  fhe  me¬ 
dium,  wifh  fhe  rafe  of  decrease  dependenf  on  fhe 
absorpfion  coefficienf.  Whaf  needs  fo  be  known 
are  fhe  absorpfion  coefficienfs  for  fhe  primary 
componenfs  of  sea  ice:  ice,  brine  and  air.  Equa¬ 
fion  3  implies  fhaf  absorpfion  coefficienfs  can  be 
defermined  by  measuring  fhe  incidenf  radiance, 
fhe  fransmiffed  radiance,  and  fhe  fhickness  of  a 
homogeneous  sample  fhaf  is  free  of  scafferers 
(Grenfell  and  Perovich  1981). 

Absorpfion  in  fhe  air  inclusions  in  sea  ice  is 
negligible,  so  absorpfion  coefficienfs  for  air  are 
assumed  fo  be  zero.  Specfral  absorpfion  coeffi¬ 
cienfs  from  fhe  ulfraviolef  fo  fhe  near-infrared 
for  ice  and  seawafer  are  shown  in  Eigure  5.  Ab¬ 
sorpfion  coefficienfs  for  ice  were  defermined 
using  pure,  bubble-free,  fresh  ice  (Grenfell  and 
Perovich  1984,  Perovich  and  Govoni  1991),  and 
absorpfion  coefficienfs  for  brine  faken  from  mea- 
suremenfs  of  clear  Arctic  wafer  (Tyler  and  Smifh 
1970,  Smifh  and  Baker  1981).  The  minimum  ab¬ 
sorpfion  and  fherefore  maximum  fransmission 
for  ice  is  in  fhe  blue  parf  of  fhe  elecfromagnefic 
specfrum  af  470  nm.  Specfral  changes  in  ab¬ 
sorpfion  coefficienf  are  exfremely  large,  span¬ 
ning  several  orders  of  magnifude  from  250  fo 
1400  nm.  Specfral  differences  in  fhe  fransmiffed 


radiance,  I  (z,  X),  are  even  more  pronounced, 
being  exponenfially  greafer  fhan  fhe  changes  in 
absorpfion  coefficienf. 

Examining  fhe  e-folding  lengfh  gives  a  better 
undersfanding  of  fhe  absorpfion  coefficienfs.  The 
e-folding  lengfh  is  fhe  amounf  of  ice  needed  fo 
reduce  fhe  incidenf  lighf  (I  (0,  A,))  by  1/e  (i.e.,  fhe 
fransmiffed  lighf  is  37%  of  fhe  incidenf).  e-folding 
lengfhs  for  ice  decrease  sharply  from  24  m  af  470 
nm,  fo  8  m  af  600  nm,  fo  2  m  af  700  nm,  fo  0.05  m 
af  1000  nm,  fo  0.006  m  af  1400  nm.  This  indicafes 
fhaf  ice  is  quife  fransparenf  in  fhe  blue,  while  if 
fakes  only  a  few  cenfimefers  of  ice  fo  absorb  mosf 
of  fhe  lighf  beyond  1000  nm. 

As  Eigure  5  indicafes,  absorpfion  coefficienfs 
for  clean  Arctic  wafer  are  similar  in  magnifude 
and  specfral  shape  fo  values  for  pure  ice.  An  ab¬ 
sorpfion  coefficienf  for  sea  ice  is  defermined 
by  combining  fhe  absorpfion  coefficienfs  for  fhe 
consfifuenf  componenfs  of  brine  and  ice  using 

fcsi  =  Vifci-rVbkb  (4) 

where  Vj  and  Vb  represenf  fhe  volume  fraction  of 
ice  and  brine,  and  k^  and  kb  ^re  fhe  absorpfion 


Wavelength  (nm) 

Figure  5.  Absorption  coefficients  of  pure,  bubble-free 
ice  (Grenfell  and  Perovich  1981,  Perovich  and  Govoni 
1991)  and  clear  sea  water  (Tyler  and  Smith  1970,  Smith 
and  Baker  1981). 
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coefficients  of  ice  and  brine  (Grenfell  1983).  Ab¬ 
sorption  by  air  is  assumed  to  be  negligible. 

Equation  4  provides  a  simple  means  of  gener¬ 
ating  an  absorption  coefficient  for  sea  ice  from 
physically  determinable  quantities  and  known 
values  of  ice  and  brine  absorption  coefficients. 
Unfortunately,  in  nature,  sea  ice  is  often  more 
than  a  combination  of  ice,  brine  and  air.  For  ex¬ 
ample,  particulates,  sediments,  ice  biota  and  dis¬ 
solved  organics  can  be  present.  If  these  impuri¬ 
ties  are  present  in  sufficient  quantity,  then  their 
absorptive  properties  must  also  be  considered.  In 
general  these  impurities  are  strongly  absorbing 
and  weakly  scattering.  Absorption  coefficients  for 
sediments  and  ice  biota  vary  depending  on  their 
composition.  Examples  of  absorption  coefficients 
for  ice  biota  (Arrigo  et  al.  1991)  are  shown  in 
Figure  6.  The  spectral  shapes  of  these  absorption 
coefficients  are  quite  different  than  those  of  ice  or 
brine.  If  sediment  or  ice  biota  are  present  in  suffi¬ 
cient  quantity,  they  should  be  explicitly  treated  in 
the  theoretical  formulation  by  modifying  eq  4. 

Scattering 

Sea  ice  is  not  a  monolithic  slab  of  pure  ice.  It 
has  an  intricate  structure  consisting  of  an  ice  ma- 


Figure  6.  Absorption  coefficients  of  biota  found  in 
congelation  ice  and  frazil  ice  (from  Arrigo  et  al.  1991). 


trix  with  inclusions  of  brine,  air  and  perhaps  solid 
salts.  Since  these  inclusions  have  different  indices 
of  refraction  than  the  surrounding  ice,  they  scat¬ 
ter  light.  The  larger  the  difference  in  index  of 
refraction  between  the  inclusion  and  the  ice,  the 
stronger  the  scattering.  Sea  ice  has  an  abundance 
of  brine  pockets  and  air  bubbles  and  therefore  is 
a  highly  scattering  medium.  In  certain  cases,  par¬ 
ticulates,  sediment,  and  ice  biota  contribute  to 
scattering,  but  air  bubbles  and  brine  pockets  are 
the  primary  scatterers  in  sea  ice  and  are  the  focus 
of  this  discussion. 

Scattering  results  from  differences  in  the  real 
indices  of  refraction  (n)  between  ice  {n  ~  1.31)  and 
the  inclusions.  With  a  greater  difference  in  index 
of  refraction,  air  bubbles  {n  ~  1.0)  are  more  strong¬ 
ly  scattering  than  brine  pockets.  The  real  part  of 
the  index  of  refraction  for  brine  depends  on  tem¬ 
perature,  increasing  from  1.34  at  -2°C  to  1.40  at 
-32°C  (Maykut  and  Light  1995).  If  the  ice  is  cold 
enough  that  solid  salts  form,  scattering  increases 
significantly,  since  these  salts  are  very  effective 
scatterers  (Perovich  and  Grenfell  1981).  The  scat¬ 
tering  coefficient  depends  not  only  on  the  amount 
of  brine  and  air,  but  on  how  it  is  distributed.  This 
complicates  matters  since  the  readily  determined 
brine  and  air  volumes  are  not  sufficient  to  define 
scattering.  The  more  difficult  to  obtain  size  distri¬ 
bution  of  the  inclusions  is  also  needed.  More  in¬ 
clusions  in  the  ice  results  in  more  scattering  and 
consequently  a  larger  scattering  coefficient.  Scat¬ 
tering  coefficients  in  sea  ice  are  large,  with  values 
typically  greater  than  10  m“^  for  warm  ice  and 
greater  than  200  m“^  for  ice  with  abundant  air 
bubbles  or  ice  colder  than  -24°G  with  precipi¬ 
tated  hydrohalite  present  (Perovich  and  Grenfell 
1982). 

Scattering  is  defined  by  two  parameters:  the 
scattering  coefficient  and  the  phase  function.  The 
scattering  coefficient  (a)  is  analogous  to  the  ab¬ 
sorption  coefficient  and  is  a  measure  of  the 
amount  of  scattering  per  unit  length.  The  phase 
function  [p(p,p',(l),  (|)')]  describes  the  angular  de¬ 
pendence  of  scattering  and  usually  is  normalized 
so  that  its  integral  over  the  full  range  of  p  and  (|)  is 
equal  to  one. 

Because  scattering  depends  on  the  intricate  and 
highly  variable  microstructure  of  sea  ice,  it  is  not 
possible  to  formulate  a  simple,  all-encompassing 
equation  to  define  the  scattering  coefficient  and 
the  phase  function,  as  we  could  for  the  absorp¬ 
tion  coefficient.  Gomplicating  matters  even  fur¬ 
ther  is  the  fact  that  in  a  highly  scattering  medium 
such  as  sea  ice,  scattering  coefficients  and  phase 


6 


functions  are  extremely  difficult  to  measure.  There 
is,  however,  one  simplifying  aspect  to  scattering 
in  sea  ice  in  the  optical  regime:  it  can  be  assumed 
to  be  independent  of  wavelength.  The  wavelength 
dependence  of  the  real  portion  of  the  index  of 
refraction  for  ice,  brine  and  air  is  very  weak  at 
optical  wavelengths  and  t5^ically  is  assumed  to 
be  constant  with  wavelength  (Grenfell  1983, 1991). 
Optical  wavelengths  are  on  the  order  of  tenths  of 
a  micrometer  to  micrometers.  The  inclusions  in 
sea  ice  have  sizes  on  the  order  of  tenths  of  a 
millimeter  for  brine  pockets  to  millimeters  for  air 
bubbles.  Since  the  scatterers  are  much  bigger  than 
the  wavelength  and  the  scatterers  are  far  apart, 
contributions  due  to  diffraction  and  interference 
can  be  ignored  (Grenfell  1983,  Bohren  and  Huff¬ 
man  1983).  The  result  of  the  weak  wavelength 
dependence  of  n,  and  the  fact  that  the  size  of  the 
scatterers  is  much  larger  than  the  wavelength,  is 
that  scattering  coefficients  and  phase  functions 
for  sea  ice  can  be  assumed  to  be  constant  with 
wavelength  (Grenfell  1983,  1991,  Perovich  1993). 
A  similar  argument  is  made  when  analyzing 


Figure  7.  Observed  (Grenfell  and  Hedrick  1983)  and 
calculated  (from  Light  1995)  phase  functions  for  sea 
ice.  0  °  is  forward  scattering  and  180°  is  backward  scat¬ 
tering:  a)  observations  of  ice  grown  at  -30  °C,  b)  obser¬ 
vations  of  ice  grown  at  -10  °C,  and  c)  calculated  esti¬ 
mates  for  ice  at  -30  °C. 


scattering  in  snow  (Bohren  and  Barkstrom  1974, 
Wiscombe  and  Warren  1980).  A  thorough  general 
discussion  of  scattering  can  be  found  in  van  de 
Hulst  (1981)  and  Bohren  and  Huffman  (1983). 

Observations  of  scattering  parameters  are  lim¬ 
ited.  Perovich  and  Grenfell  (1982)  estimated  scat¬ 
tering  coefficients  for  young  ice  from  observa¬ 
tions  of  albedo  and  transmittance.  They  found 
that  scattering  coefficients  ranged  from  8.9  m“^ 
for  melting  young  ice  to  19.6  m“^  for  cold  young 
ice  to  420  m“^  for  very  cold  young  ice  with  pre¬ 
cipitated  solid  salts  present.  Grenfell  and  Hedrick 
(1983)  used  small  samples  of  young  ice  to  mea¬ 
sure  phase  functions  for  sea  ice.  Phase  functions 
for  columnar  ice  samples  grown  at  -10°G  and 
-30°G  are  shown  in  Figure  7.  The  phase  function 
is  strongly  forward-peaked,  with  forward  scat¬ 
tering  being  more  than  a  factor  of  50  greater  than 
side  or  backward  scattering.  However,  although 
small  samples  were  used,  there  was  still  multiple 
scattering,  and  consequently  the  results  represent 
only  an  approximation  to  the  true  single  scatter¬ 
ing  albedo  and  phase  function.  Multiple  scatter¬ 
ing  tends  to  smooth  and  reduce  the  angular  de¬ 
pendence  of  the  measured  phase  function. 

Numerical  calculations  have  been  used  to 
supplement  the  relatively  sparse  observational 
data  (Grenfell  1983,  1991).  Phase  functions  are 
calculated  using  a  Mie  scattering  model  with  the 
indices  of  refraction  for  ice  and  brine  and  inclu¬ 
sion  size  distributions  as  input  parameters  (Bohren 
and  Huffman  1983).  A  calculated  phase  function 
for  sea  ice  at  -30  °G  with  brine  pockets  with  a 
radius  of  0.02  mm  (Light  1995)  is  compared  to 
observed  values  in  Figure  7.  As  expected,  the  cal¬ 
culated  phase  function  is  more  strongly  forward- 
peaked  than  the  multiply-scattered  observed. 

Though  we  do  not  have  a  quantitative  under¬ 
standing  of  the  relationship  between  scattering 
and  ice  physical  properties,  a  qualitative  grasp  is 
sufficient  for  our  purposes.  To  interpret  observa¬ 
tions  of  optical  properties  the  important  theoreti¬ 
cal  points  are  1)  absorption  coefficients  for  ice 
and  brine  depend  strongly  on  wavelength,  2)  scat¬ 
tering  coefficients  and  phase  functions  for  sea  ice 
are  constant  with  wavelength,  3)  increasing  the 
number  of  inclusions  in  sea  ice  increases  the 
amount  of  scattering,  4)  air  bubbles  scatter  more 
strongly  than  brine  pockets,  and  5)  scattering  in 
sea  ice  is  strongly  forward  peaked.  With  this 
theoretical  foundation  regarding  the  underly¬ 
ing  physics  of  radiative  transfer  in  sea  ice,  it  is 
time  to  revisit  the  question  of  "What  is  the  albedo 
of  sea  ice?" 
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OBSERVATIONS 

There  is  a  large  observational  dataset  of  sea  ice 
optical  properties,  particularly  of  sea  ice  albedos 
(Perovich  ef  al.  1986).  In  fhis  secfion  we  presenf 
an  overview  of  fhese  observafions  in  fhe  confexf 
of  illusfrafing  how  fhe  opfical  properfies  of  sea 
ice  are  affecfed  by  fhe  physical  properfies.  We 
invesfigafe  fhe  effecfs  of  ice  type,  surface  condi¬ 
tions,  ice  fhickness,  ice  brine  volume,  and  impu¬ 
rities  on  albedo,  reflecfance,  fransmiffance,  and 
extinction  coefficienf.  The  simplifying  beaufy  of 
opfical  properly  observafions  is  fhaf,  af  leasf  from 
400  nm  fo  750  nm,  whaf  you  see  is  whaf  you  gef. 
If  fhe  ice  looks  while,  fhen  ifs  albedo  will  be  high 
and  relafively  consfanf  wifh  wavelengfh.  Simi¬ 
larly,  fhe  specfral  albedo  of  a  blue-looking  melf 
pond  will  have  a  peak  befween  400  and  500  nm. 

Albedos 

Albedos  are  sensitive  fo  fhickness  during  fhe 
initial  sfages  of  ice  growfh.  Weller  (1972)  mea¬ 
sured  fofal  albedos  in  a  freezing  lead  and  found  a 
rapid  rise  in  albedo  from  0.08  fo  0.40  as  fhe  ice 
grew  from  open  wafer  fo  a  fhickness  of  0.30  m, 
followed  by  a  more  gradual,  asympfofic  increase 
as  fhe  ice  confinued  fo  grow.  For  ice  fhicker  fhan 
approximafely  0.8  m  fofal,  albedo  shows  little 
change  wifh  fhickness  (Maykuf  1982). 

Specfral  changes  in  albedo  during  inifial  ice 
growfh  are  ploffed  in  Figure  8.  The  ice  was  grown 
in  fhe  laboratory  af  a  consfanf  air  femperafure  of 
-20°C  and  had  a  columnar  crysfal  sfrucfure 
(Perovich  1979).  Albedo  increased  wifh  fhickness 
af  all  wavelengfhs.  As  fhe  ice  grows  fhicker,  op- 
porfunifies  for  backscaffering  in  fhe  ice  are  added 
and  af  firsf  albedos  rise  rapidly.  However,  fhe 
pafh  lengfh  of  fhe  backscaffered  lighf  also  in¬ 
creases,  until  finally  only  a  negligible  amounf  of 
lighf  penefrafes  fo  fhe  boffom  of  fhe  ice,  is 
backscaffered,  and  emerges  from  fhe  ice  surface 
wifhouf  being  absorbed.  Af  fhis  poinf  increasing 
fhe  ice  fhickness  no  longer  direcfly  affecfs  fhe 
albedo  and  fhe  ice  is  opfically  fhick.  As  Figure  8 
indicafes,  fhis  asympfofic  ice  fhickness  is  smaller 
af  longer  wavelengfhs.  Af  shorter  wavelengfhs 
albedos  are  still  increasing  when  fhe  ice  is  0.25  m 
fhick,  while  af  longer  wavelengfhs  (beyond  700 
nm),  fhe  asympfofic  nafure  of  fhe  albedo  increase 
is  evidenf.  This  is  a  direcf  resulf  of  fhe  increase  in 
absorption  as  wavelengfh  increases.  This  is  consis- 
fenf  wifh  our  earlier  commenf:  specfral  variations 
in  opfical  properfies  are  due  fo  absorption.  A  closer 
look  af  Figure  8  shows  fhaf  fhe  rapid  rise  in  al- 


Figure  8.  Laboratory  observations  of  the  increase  in 
spectral  albedo  during  initial  ice  growth  (Perovich 
1979).  The  ice  was  grown  at  an  air  temperature  of 
-20  °C. 

bedo  did  nof  begin  until  fhe  ice  was  0.05  m  fhick. 
During  fhe  firsf  0.05  m  of  growfh,  fhe  ice  in  fhis 
experimenf  had  nof  yef  begun  fo  cool  and  brine 
volumes  were  quite  large.  Because  of  fhis  fhere 
was  liffle  scaffering  in  fhe  ice.  As  fhe  ice  cooled 
and  fhe  large  brine  pockefs  fragmented  info  many 
smaller  pockefs,  fhere  was  more  scaffering  and 
fhe  albedo  increased.  This  observafion  illusfrafes 
again  fhaf  fhe  opfical  properfies  of  sea  ice  are 
often  more  complicated  fhan  we  would  expecf. 

Of  course,  even  for  fhick  ice,  sea  ice  albedos 
vary.  From  fhe  previous  discussion  we  can  see 
fhaf  albedo  is  sensitive  fo  fhe  ice  surface  condi- 
fions.  Specfral  albedos  for  mulfiyear  ice  are  plof- 
fed  in  Figure  9  (faken  from  Grenfell  and  Maykuf 
1977).  These  albedos  represenf  a  possible  evolu¬ 
tionary  sequence  from  spring  fo  summer  as  melf 
occurs,  and  fhe  ice  cover  changes  from  snow- 
covered  ice  fo  bare  ice  or  frozen  ponds  fo  melfing 
ice  fo  ponded  ice.  Snow  albedos  (curve  a)  are 
large  (-0.9)  and  nearly  consfanf  wifh  wavelengfh 
in  fhe  visible;  fhe  snow  appears  brighf  and  whife. 
Scaffering  coefficienfs  for  snow  are  so  large  fhaf, 
in  fhe  visible,  absorption  has  liffle  impacf  on  fhe 
albedo  and  fhere  is  no  wavelengfh  dependence. 
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A  0.1-m-thick  layer  of  wind-packed  Arctic  snow 
is  sufficient  to  eliminate  any  contribution  to  the 
albedo  from  the  underlying  ice. 

Spectral  albedos  decrease  as  the  surface 
changes  from  snow-covered  ice  to  cold,  bare 
multiyear  ice  (curve  b).  White  multiyear  ice  typi- 


Wavelength  (nm) 


Figure  9.  Spectral  albedos  for  a  possible  evolu¬ 
tionary  sequence  of  multiyear  ice  (Grenfell  and 
Maykut  1977):  a)  snow-covered  ice,  b)  cold  bare 
ice,  c)  melting  bare  ice,  d)  early-season  melt 
pond  and  e)  mature  melt  pond. 


cally  has  a  drained  bubbly  surface  layer  with 
plenty  of  air  bubbles  which,  while  not  as  strongly 
scattering  as  snow,  still  contributes  considerable 
scattering.  The  result  is  an  overall  decrease  in 
albedo  of  approximately  10%  and  a  slight  wave¬ 
length  dependence.  As  the  ice  warms  and  begins 
to  melt  (curve  c),  albedos  continue  to  decrease 
with  a  more  evident  wavelength  dependence.  This 
is  due  to  a  decrease  in  scattering  as  the  ice  melts 
and  some  of  the  air  voids  fill  with  water.  In  some 
areas  of  the  ice  cover,  water  collects  on  the  sur¬ 
face,  forming  melt  ponds  (curve  d).  As  the  melt 
season  progresses  these  ponds  can  get  deeper 
(curve  e).  Albedos  of  ponded  ice  are  character¬ 
ized  by  a  maximum  in  the  400-500  nm  region 
and  a  precipitous  decrease  between  500  and  800. 
The  melt  ponds  look  blue.  This  spectral  behavior 
is  due  to  the  transparency  of  the  water  at  shorter 
wavelengths;  albedos  below  500  nm  are  deter¬ 
mined  primarily  by  the  scattering  properties  of 
the  underlying  ice.  From  500  to  800  nm  the  al¬ 
bedo  becomes  increasingly  insensitive  to  the  un¬ 
derlying  ice  as  the  absorption  in  the  water  be¬ 
comes  the  dominant  factor.  Above  800  nm 
absorption  in  the  water  is  so  great  that  pond  albe¬ 
dos  are  essentially  determined  by  Fresnel  reflec¬ 
tion  at  the  surface  and  are  independent  of  wave¬ 
length.  These  results  indicate  that  both  the 
magnitude  and  the  shape  of  the  spectral  albedos 
are  extremely  sensitive  to  the  amount  of  liquid 
water  present  in  the  upper  part  of  the  ice. 

When  skies  are  clear,  the  wavelength  interval 
from  1000  to  2500  nm  can  contain  up  to  25%  of  the 
total  incident  shortwave  energy  (Grenfell  and 
Perovich  1984),  so  albedos  in  this  re¬ 
gion  can  have  a  significant  impact  on 
the  heat  and  mass  balance  of  the  ice. 
Figure  10  shows  a  spectral  albedo  se¬ 
quence  that  first-year  ice  might  follow 
as  it  progresses  through  a  melt  cycle 
from  (a)  ice  covered  by  cold  dry  snow, 
(b)  to  ice  covered  by  melting  snow,  (c)  to 
bare  ice  with  a  crumbly  surface  layer, 
and  (d)  to  melting  first-year  blue  ice. 
Concentrating  on  wavelengths  beyond 
1000  nm,  a  continual  downward  trend 
is  evident,  with  albedo  reaching  a  mini¬ 
mum  at  about  1500  nm.  Local  maxima 
are  located  at  1350,  1900  and  2300  nm 
and  correspond  to  minima  in  the  ab¬ 
sorption  spectrum  for  ice.  In  general, 
sea  ice  and  snow  albedos  at  longer 
wavelengths  are  significantly  smaller 
than  values  at  visible  wavelengths. 


Figure  10.  A  spectral  albedo  sequence  that  first-year  ice  might 
follow  through  a  melt  cycle  (Grenfell  and  Perovich  1984):  a)  ice  cov¬ 
ered  by  cold  dry  snow,  to  b)  ice  covered  by  melting  snow,  to  c)  bare 
ice  with  a  crumbly  surface  layer,  to  (d)  melting  first-year  blue  ice. 
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Comparing  albedos  in  Figures  9  and  10  dem¬ 
onstrates  that,  for  equivalent  conditions,  multiyear 
ice  albedos  are  typically  larger  than  first-year  ice 
values.  Multiyear  ice  has  undergone  a  summer 
melt  season,  with  the  attendant  surface  melting 
and  brine  drainage.  This  resulfs  in  a  well-devel¬ 
oped  surface-scaffering  layer  wifh  many  air 
bubbles. 

As  a  resulf  of  fhe  decrease  in  albedo  af  longer 
wavelengfhs,  fofal  albedos  (aj)  are  greafer  under 
cloudy  skies  fhan  under  clear  skies.  The  fofal  al¬ 
bedo  depends  on  fhe  specfral  albedo  and  fhe  spec- 
fral  incidenf  irradiance  (eq  1).  Clouds  absorb  more 
sfrongly  in  fhe  infrared  fhan  in  fhe  visible.  There¬ 
fore  on  cloudy  days  a  greafer  porfion  of  fhe  inci¬ 
denf  irradiance  is  af  visible  wavelengfhs,  where 
fhe  albedo  is  larger.  Tofal  albedos  under  cloudy 
skies  are  f5q)ically  8-12%  larger  fhan  clear  sky 
values  (Grenfell  and  Maykuf  1977,  Grenfell  and 
Perovich  1984).  When  fhe  incidenf  direcf  beam 
componenf  is  significanf,  bofh  specfral  and  fofal 
albedos  increase  as  fhe  solar  zenifh  angle  increases 
(sun  closer  fo  fhe  horizon)  due  fo  enhanced  specu¬ 
lar  reflection  (Perovich  and  Grenfell  1982)  and  fo 
forward  scaffering  allowing  fhe  phofons  fo  es¬ 
cape  fhe  medium  fasfer. 

Surface  conditions  have  a  sfrong  impacf  on 
albedo,  buf  fhe  infernal  sfafe  and  sfrucfure  of  fhe 
ice  are  also  significanf.  As  we  have  seen,  fhe  pres¬ 
ence  of  air  bubbles  in  fhe  upper  porfion  of  fhe  ice 
enhances  albedo.  Brine  volume  is  anofher  impor- 
fanf  ice  physical  properly  fhaf  we  mighf  expecf  fo 
have  some  impacf  on  albedo.  Perovich  and 
Grenfell  (1981)  invesfigafed  fhe  influence  of  brine 
volume  on  albedo  for  young  ice.  Resulfs  from 
fhree  laboratory  experimenfs  are  summarized  in 
Figure  11.  In  each  experimenf  fhe  ice  was  grown 
af  a  selected  air  femperafure  (-37°,  -30°  and 
-10°G),  and  fherefore  a  differenf  growfh  rate,  fo  a 
fhickness  of  approximafely  0.25  m.  Brine  volume 
was  fhen  varied  by  warming  fhe  ice.  For  each 
experimenf  fhere  is  a  continual  decrease  in  al¬ 
bedo  as  brine  volume  increases.  As  fhe  ice  warmed 
and  fhe  brine  volume  increased,  individual  brine 
pockefs  coalesced,  forming  larger  buf  fewer  in¬ 
clusions.  Thus,  fhe  resulf  of  fhe  warming  was  a 
reducfion  in  fhe  number  of  brine  pockefs  and  in 
fhe  amounf  of  scaffering.  Gomparing  befween  ex¬ 
perimenfs,  we  also  see  fhaf  af  a  given  brine  vol¬ 
ume  fhere  is  considerable  variabilify  in  fhe  ob¬ 
served  albedo,  wifh  fasfer  grown  ice  (tower  air 
femperafure)  having  larger  albedos  (Perovich  and 
Grenfell  1981).  More  rapidly  grown  ice  has  smaller 
plafelef  and  crysfal  sizes  and  more  brine  inclu- 


Brine  Volume  (%) 


Figure  11.  Observations  of  total  albedo  vs.  brine  volume 
for  young  ice  (from  Perovich  and  Grenfell  1981).  Three 
experiments  were  performed  where  the  ice  was  grown  to 
approximately  0.25  m  thick  and  then  warmed,  in  stages, 
to-2°C. 


sions  (Weeks  and  Plamilfon  1962,  Lofgren  and 
Weeks  1969,  Weeks  and  Ackley  1982).  This  leads 
fo  fhe  imporfanf  conclusion  fhaf  nof  only  is  fhe 
volume  of  brine  imporfanf,  buf  how  if  is  disfrib- 
ufed  is  also  significanf.  For  a  given  volume  of 
brine,  fhere  is  more  scaffering  if  fhaf  brine  is  dis- 
fribufed  info  many  small  brine  inclusions,  rafher 
fhan  a  single  large  one.  The  same  conclusion  is 
frue  for  air  bubbles. 

All  of  fhe  albedos  presenfed  so  far  have  been 
for  Arcfic  sea  ice.  Are  albedos  for  Anfarcfic  ice 
differenf?  Specfral  albedos  for  young  sea  ice 
grown  off  of  Easf  Anfarcfica  are  ploffed  in  Figure 
12  (Allison  ef  al.  1993).  These  albedos  show  fhe 
same  general  properties  as  Arcfic  sea  ice  resulfs; 
an  increase  as  fhe  ice  grows  fhicker  and  a  gradual 
wavelengfh  dependence  wifh  larger  albedos  af 
shorfer  wavelengfhs  (Schlosser  1988,  Allison  ef 
al.  1993).  There  are  differences  in  ice  sfrucfure 
befween  Anfarcfic  and  Arcfic  sea  ice.  Anfarcfic 
sea  ice  has  much  more  frazil  ice  fhan  Arcfic  sea 
ice  and  is  somewhaf  more  saline  (Gow  and  Tucker 
1990).  Surface  condifions  also  differ  wifh  signifi¬ 
canf  amounfs  of  flooded  snow-covered  ice  buf 
very  little  ponded  ice  in  fhe  Anfarcfic  (Andreas 
and  Ackley  1982).  Because  of  fhis  fhere  may  be 
differences  in  optical  properties  befween  Anfarc¬ 
fic  and  Arcfic  sea  ice,  buf  any  differences  between 
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Reflectance  Factor 


the  two  cases  are  smaller  than  differences  within 
fhe  fwo. 

Reflectance 

When  considering  light  reflected  from  sea  ice, 
fhe  albedo  is  fhe  paramefer  of  prime  climatologi¬ 
cal  imporfance.  However,  opfical  remofe  sensing 


Figure  12.  Spectral  albedos  of  Antarctic  sea  ice  (from 
Allison  et  al.  1993). 


insfrumenfs  mounfed  on  aircraff  or  safellifes  typi¬ 
cally  have  narrow  fields  of  view  and  measure 
reflecfed  radiance,  rafher  fhan  irradiance.  To  con- 
verf  fhese  observafions  of  radiance  to  an  esfimafe 
of  irradiance,  informafion  on  fhe  bidirecfional  re- 
flecfance  disfribufion  funcfion  (BRDF)  is  needed. 
Under  cloudy  conditions  fhe  incidenf  radiation 
field  is  diffuse,  so  fhe  lighf  reflecfed  from  fhe 
surface  is  also  diffuse.  Reflecfed  radiance  is  es¬ 
sentially  fhe  same  in  any  direcfion  and  fhe  re¬ 
flecfed  irradiance  is  easily  compufed.  When  if  is 
sunny,  fhough,  fhe  incidenf  radiafion  consisfs  of 
a  diffuse  sky  componenf  plus  a  very  sfrong  solar 
direcf  beam,  and  fhe  incidenf  radiance  field  is 
sfrongly  anisofropic.  For  fhese  incidenf  condifions 
fhe  angular  disfribufion  of  reflecfed  radiance  can 
be  complex.  Angular  reflecfances  (reflecfed  radi¬ 
ance  normalized  fo  a  white  reference  sfandard) 
measured  under  surmy  skies  for  snow-covered 
ice  and  bare  blue  ice  are  plotted  in  Figure  13. 
Snow-covered  ice  reflecfances  are  fairly  consfanf 
wifh  angle,  excepf  for  a  30%  increase  af  fhe  angle 
of  reflecfion  of  fhe  solar  beam.  The  peak  in  reflec- 
fance  af  fhe  angle  of  reflecfion  is  even  more  pro¬ 
nounced  in  fhe  bare  blue  ice  case,  wifh  an  in¬ 
crease  fo  nearly  fwice  fhe  value  of  fhe  albedo.  Af 
ofher  angles,  blue  ice  reflecfances  are  equal  fo  or 
slighfly  less  fhan  fhe  albedo.  The  differences  be- 
fween  R  and  a  show  fhe  imporfance  of  fhe  BRDF 
in  defermining  albedos  from  observafions  of  re¬ 
flecfed  radiance.  The  presence  of  any  sysfemafic 
fopographic  feafures,  such  as  sasfrugi,  will  fur- 
fher  complicate  fhe  BRDF. 


Figure  13.  Bidirectional 
reflectance  distribution 
function  at  450  nm  for 
snow-covered  ice  and  bare 
blue  ice  (from  Perovich 
1994).  Rq  is  the  normal¬ 
ized  reflected  radiance  at 
nadir  and  a  is  the  albedo. 
The  measurements  were 
made  at  a  solar  zenith 
angle  of  60°  under  clear 
skies. 
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The  peak  in  reflectance  is  largely  due  to  specu¬ 
lar  reflection  of  the  direct  solar  beam.  Specular 
reflection  is  light  reflected  from  the  surface  of  the 
medium  in  the  direction  of  the  angle  of  reflec¬ 
tion.  Its  magnitude  depends  on  the  angle  of  inci¬ 
dence  and  the  index  of  refraction  according  to: 


4(e)  =  ^o(e)f 


sin^  (9-(|)’^) 
sin^  (0-i-(|)’^) 


-I- 


tan^ 

tan^  (6-i-(|)’^)^ 


where  Iq{Q)  =  the  incident  direct  solar  beam  ra¬ 
diance, 

tj.(0)  =  the  reflected  radiance 

0  =  the  zenith  angle  of  incident  and 
reflected  radiance 
0*=  [n  arcsin(0)]“i 

n  =  the  index  of  refraction  of  the  me¬ 
dium  (Bom  and  Wolfe  1965). 

As  the  zenith  angle  increases  (sun  gets  closer  to 
the  horizon),  the  specular  reflection  increases.  The 
specular  component  is  larger  for  smooth  surfaces 
such  as  blue  ice  or  melt  ponds  and  smaller  for 
rough  surfaces  such  as  snow  or  drained  white 
ice. 


Transmission 

The  magnitude  and  spectral  distribution  of 
light  transmitted  through  the  ice  cover  depends 
on  the  physical  composition  of  the  ice,  the  thick¬ 
ness  of  the  ice,  and  the  surface  conditions.  As 
thickness  increases,  light  transmission  through  the 
ice  drops  off  roughly  exponentially.  The  influ¬ 
ence  of  surface  conditions  on  light  transmission 
through  the  ice  is  illustrated  in  Figure  14  (Maykut 
and  Grenfell  1975).  Even  a  thin  (0.25-m)  layer  of 
highly  scattering  snow  can  reduce  transmittances 
through  the  ice  cover  to  less  than  1%  (curve  a).  As 
the  snow  melts,  scattering  decreases  and  the  trans¬ 
mittance  increases  (curve  b)  until  the  snow  is  gone 
(curve  c).  The  presence  of  melt  ponds  greatly  re¬ 
duces  scattering  in  the  upper  portion  of  the  ice 
and  enhances  transmission  in  the  visible  (curve 
d).  Light  levels  beneath  ponded  first-year  ice  are 
at  least  a  factor  of  three  greater  than  those  be¬ 
neath  white  ice  of  the  same  thickness  (Grenfell 
and  Maykut  1977). 

Sea  ice  can  be  a  prime  habitat  for  biological 
organisms.  The  growth  of  these  organisms  is  in¬ 
fluenced  by  the  amount  of  light  available  in  and 
under  the  ice  (Soo  Hoo  et  al.  1987,  Gota  and  Horne 
1989,  Arrigo  et  al.  1991).  Just  as  the  ice  biota  are 
affected  by  the  light  levels,  they  in  turn  can  re¬ 


Figure  14.  The  influence  of  surface  conditions 
on  light  transmission  (from  Maykut  and  Grenfell 
1977).  In  all  cases  the  ice  thickness  was  1.85  m. 
Surface  conditions  were  a)  blue  ice  covered  by 
0.25  m  of  melting  snow,  b)  blue  ice  covered  by 
0.12  m  of  melting  snow,  c)  white  ice  and  d)  blue 
ice  covered  by  a  0.05-m  melt  pond. 

duce  light  transmission  and  change  its  spectral 
composition  (Maykut  and  Grenfell  1975,  Soo  Hoo 
et  al.  1987,  Perovich  et  al.  1993).  As  an  example  of 
this,  spectral  transmittances  for  1.5-m-thick  first- 
year  ice  with  a  0.05-m  snow  cover  and  with  a 
0.19-m  snow  cover  are  plotted  in  Figure  15.  Sur¬ 
prisingly,  the  transmission  is  less  under  the  thin¬ 
ner  snow  cover.  This  is  a  direct  result  of  a  50% 
higher  algal  biomass  (157  vs.  117  mg  chlorophyll 
m“^)  at  the  thin  snow  site.  In  addition  to  the  over¬ 
all  reduction,  the  presence  of  the  additional  bio¬ 
mass  results  in  enhanced  losses  in  the  blue  end  of 
the  spectrum  and  a  pronounced  drop  in  transmit¬ 
tance  at  670  nm. 

Extinction  coefficient 

A  more  fundamental  way  to  characterize  light 
penetration  through  sea  ice  is  by  an  extinction 
coefficient  (k;)^^).  The  extinction  coefficient  is  a  mea¬ 
sure  of  loss  due  to  scattering  and  absorption  and 
is  typically  determined  from  measurements  of 
incident,  reflected  and  transmitted  light.  Trans¬ 
mission  measurements  are  difficult  and  demand¬ 
ing,  and  consequently  there  have  been  far  fewer 


12 


Figure  15.  Observed  spectral  transmittances  for  1.5- 
m-thick first-year  ice  with  a)  0.05-m  snow  cover  plus 
157  mg  chlorophyll  m~^  biomass  and  b)  0.19  m  snow 
cover  plus  117  mg  chlorophyll  mr^. 

observations  of  extinction  coefficient  than  of  al¬ 
bedo.  Most  of  the  reported  sea  ice  spectral  extinc¬ 
tion  coefficients  have  been  calculated  using  a  two- 
stream  radiative  transfer  model  (Grenfell  and 
Maykut  1977,  Perovich  and  Grenfell  1981). 

Figure  16  summarizes  spectral  extinction  coef¬ 
ficients  culled  from  a  number  of  sources  for  nine 
distinct  cases:  dry  snow,  melting  snow,  ice  below 
the  eutectic  point  with  solid  salts  present,  the 
surface  scattering  layer  of  white  ice,  the  interior 
of  white  ice,  cold  blue  ice,  melting  blue  ice,  bubble- 
free  fresh  ice,  and  clear  Arctic  water  (Grenfell  and 
Maykut  1977,  Perovich  and  Grenfell  1981,  Smith 
and  Baker  1981).  The  range  of  over  one  to  two 
orders  of  magnitude  in  the  extinction  coefficients 
shows  the  tremendous  variation  in  attenuation 
between  different  snow  and  ice  types.  Sea  ice  and 
snow  curves  all  show  relatively  constant  values 
in  the  400-  to  500-nm  region,  followed  by  strongly 
increased  attenuation  at  longer  wavelengths. 
Again,  as  was  the  case  for  albedo,  the  magnitude 
of  the  extinction  coefficient  is  largely  a  function 
of  the  amount  of  scattering,  while  the  wavelength 
dependence  is  determined  by  absorption.  The 
greatest  attenuation  occurs  in  cold  snow,  where 


Figure  16.  Spectral  extinction  coefficients  for  nine 
distinct  cases:  a)  dry  snow  (Grenfell  and  Maykut 
1977),  b)  ice  below  the  eutectic  point  with  solid  salts 
present  (Perovich  and  Grenfell  1981),  c)  melting  snow 
(Grenfell  and  Maykut  1977),  d)  surface  scattering 
layer  of  white  ice  (Grenfell  and  Maykut  1977),  e)  the 
interior  of  white  ice  (Grenfell  and  Maykut  1977),  f) 
cold  blue  ice  (Grenfell  and  Maykut  1977),  g)  melting 
blue  ice  (Grenfell  and  Maykut  1977),  h)  bubble-free 
fresh  ice  (Grenfell  and  Perovich  1981),  and  i)  clear 
Arctic  water  (Smith  and  Baker  1981). 

spectral  extinction  coefficients  are  about  20  times 
larger  than  those  of  melting  blue  ice. 

Goefficients  in  melting  snow  are  above  half 
those  in  cold  dry  snow.  Extinction  coefficients  for 
very  cold  ice  below  the  eutectic  point  are  quite 
large,  comparable  to  values  for  snow.  In  this  case 
solid  salts  precipitate  in  the  interior  of  the  sea  ice. 
These  precipitated  salts  are  small  and  plentiful 
and,  with  an  index  of  refraction  of  approximately 
1.5,  are  effective  scatterers.  The  drained  surface 
layer  of  multiyear  ice  (white  ice  scattering)  con¬ 
tains  an  abundance  of  air  inclusions  which  formed 
as  a  result  of  brine  drainage.  These  air  inclusions 
cause  considerable  scattering,  and  extinction  co¬ 
efficients  are  large.  In  the  interior  of  white  ice 
there  are  fewer  air  bubbles,  and  extinction  coeffi¬ 
cients  are  correspondingly  smaller.  In  the  blue  ice 
cases  the  inclusions  are  primarily  brine  pockets. 
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rather  than  air  bubbles,  and  the  amount  of  scat¬ 
tering  is  less  and  extinction  coefficients  are  re¬ 
duced.  The  much  smaller  values  of  exfincfion  co- 
efficienf  for  bubble-free  ice  and  clear  Arctic  wafer 
illusfrafe  how  significanf  scaffering  is  in  sea  ice. 
The  imporfance  of  scaffering  is  illusfrafed  by  fhe 
rough  rule  of  fhumb  fhaf  exfincfion  fhrough  1  cm 
of  snow  is  approximafely  fhe  same  as  fhrough  10 
cm  of  ice  or  100  cm  of  wafer. 

As  was  fhe  case  for  albedo  (Fig.  11),  exfincfion 
coefficienfs  also  depend  on  fhe  infernal  sfrucfure 
of  fhe  ice  (Zaneveld  1966,  Grenfell  and  Maykuf 
1977,  Perovich  and  Grenfell  1981,  Gilberf  and 
Bunfzen  1986).  Exfincfion  coefficienfs  decrease 
during  warming  as  fhe  brine  volume  increases 
and  fhe  number  of  inclusions  decrease.  Also,  af  a 
given  brine  volume  exfincfion  coefficienfs  are 
larger  for  fasfer  grown  ice,  which  has  more  inclu¬ 
sions.  In  fhese  experimenfs  fhe  ice  was  changing 
infernally,  buf  fhe  only  change  in  surface  condi¬ 
tions  was  a  slighf  wetting  as  fhe  air  femperafure 
approached  0°G.  The  resulfs  would  be  quife  dif- 
ferenf  if  fhere  were  brine  drainage  from  fhe  sur¬ 
face  layer  of  fhe  ice  as  a  resulf  of  fhe  warming.  In 
fhaf  case  fhe  resulfing  air  voids  would  form  a 
highly  scaffering  surface  layer,  and  albedos  and 
exfincfion  coefficienfs  would  increase.  This  would 
be  expecfed  in  fhicker  ice  wifh  more  freeboard. 
Such  an  effecf  has  been  observed  in  fhe  Anfarcfic, 
where  low  humidifies  keep  fhe  ice  surface  free  of 
wafer  during  melf  (Andreas  and  Ackley  1981). 
Observations  made  in  McMurdo  Sound,  Anfarc- 
fica  (Trodahl  ef  al.  1987,  Buckley  and  Trodahl  1987, 
Trodahl  and  Buckley  1990),  have  shown  fhaf  as 
fhe  ice  warms,  a  drained  surface  layer  forms  re¬ 
sulfing  in  an  increase  in  backscaffer  and  a  de¬ 
crease  in  fransmiffance. 

Observafions  of  fofal  lighf  fransmission  have 
been  used  fo  defermine  wavelengfh-infegrafed, 
or  fofal,  exfincfion  coefficienfs  (Kj).  Values  for  sea 
ice  are  in  fhe  1.1  fo  1.5  m“^  range  (Unfersfeiner 
1961,  Ghemogovskiy  1963,  Thomas  1963,  Weller 
and  Schwerdffeger  1967).  Exfincfion  coefficienfs 
for  snow  are  much  larger,  varying  from  4.3  m“^ 
for  dense  Anfarcfic  snow  (Weller  and  Schwerdf¬ 
feger  1967)  fo  as  high  as  40  m“^  in  freshly  fallen 
snow  (Thomas  1963).  Though  fofal  exfincfion  co¬ 
efficienfs  are  simpler  fo  measure  and  simpler  fo 
use  compufafionally  fhan  specfral  values,  fhey 
are  severely  limifed.  The  fofal  exfincfion  coeffi- 
cienf  combines  confribufions  from  differenf  wave- 
lengfhs  and  fherefore  depends  on  fhe  specfral  dis- 
fribufion  of  fransmiffed  irradiance,  which  in  furn 
depends  on  fhe  specfral  incidenf  irradiance,  fhe 


Table  1.  Values  of  ig  and  iq 
(Grenfell  and  Maykut  1977). 


Case 

Kf.  (m  h 

Clear 

Blue  ice 

0.43 

1.5 

White  ice 

0.18 

1.6 

Cloudy 

Blue  ice 

0.63 

1.4 

White  ice 

0.35 

1.5 

specfral  albedo  and  fhe  specfral  exfincfion  coeffi- 
cienf .  Since  all  of  fhese  quanfifies  vary  wifh  wave- 
lengfh,  fhe  fofal  exfincfion  coefficienf  does  nof 
depend  enfirely  on  fhe  properfies  of  fhe  ice.  As 
was  fhe  case  for  fofal  albedo,  fhe  fofal  exfincfion 
coefficienf  depends  on  sky  condifions.  On  sunny 
days  fhe  incidenf  irradiance  has  a  larger  longwave 
componenf,  which  is  absorbed  rapidly  in  fhe  ice, 
resulfing  in  higher  values  of  Kj.  More  significanfly, 
Kj  exhibifs  a  sfrong  depfh  dependence  near  fhe 
surface.  Observafions  have  shown  fhaf  fhe  spec¬ 
fral  fransmiffance  changes  greafly  near  fhe  sur¬ 
face  of  fhe  ice  due  fo  fhe  rapid  exfincfion  of  fhe 
longer  wavelengfhs.  Gorrespondingly  Kj  is  large 
near  fhe  ice  surface  and  decreases  by  more  fhan 
an  order  of  magnifude  in  fhe  fop  0.1  m  of  fhe  ice 
(Grenfell  and  Maykuf  1977).  Below  0.1  m,  only 
visible  lighf  remains,  where  fhe  specfral  depen¬ 
dence  of  is  weaker,  and  changes  in  Kj  wifh 
depfh  are  small.  Tofal  exfincfion  coefficienfs  have 
been  used  in  sea  ice  fhermodynamic  models 
(Maykuf  and  Unfersfeiner  1971)  fo  calculafe  fhe 
surface  heaf  balance  and  solar  heating  in  fhe  ice 
inferior.  To  do  fhis  Maykuf  and  Unfersfeiner  (1971) 
modified  fhe  exponenfial  decay  law  fo  fhe  form: 

(z)  =  ig  (0)e~'^‘^  for  z  >  0.1  m 

where  ig  is  fhe  fracfion  of  fhe  wavelengfh-infe¬ 
grafed  incidenf  irradiance  fransmiffed  fhrough  fhe 
fop  0.1  m  of  fhe  ice  and  k,.  is  fhe  fofal  exfincfion 
coefficienf  in  fhe  ice  below  0.1  m.  Values  of 
below  0.1  m  and  ig  defermined  from  field  obser¬ 
vafions  (Grenfell  and  Maykuf  1977)  are  summa¬ 
rized  in  Table  1.  There  is  more  scaffering  in  whife 
ice  fhan  blue  ice,  resulfing  in  a  smaller  ig  and  a 
larger  Kj. 

Beam  spread 

While  much  of  fhe  observafional  emphasis  has 
been  on  measuremenfs  of  fransmiffed  solar  irra¬ 
diance  fo  defermine  fransmiffance  and  exfincfion 
coefficienf,  measuremenfs  using  arfificial  lighf 


14 


sources  have  also  been  made.  In  particular,  stud¬ 
ies  were  conducted  examining  the  spreading  of  a 
collimated  beam  as  it  passes  through  sea  ice 
(Trodahl  et  al.  1987,  Gilbert  and  Schoonmaker 
1990,  Voss  and  Schoonmaker  1992,  Voss  et  al. 
1992).  In  these  experiments  a  collimated  beam  of 
lighf  was  incidenf  on  eifher  fhe  surface  or  boffom 
of  fhe  ice  and  fhe  spafial  disfribufion  of  fhe  emer- 
genf  irradiance  was  measured.  Examining  fhe 
peak  magnifude  and  fhe  spafial  disfribufion  of 
irradiance  provide  information  on  scaffering  and 
absorption  in  fhe  ice.  Laborafory  sfudies  indicafe 
fhaf  scaffering  in  fhe  ice  is  quife  sfrong,  wifh  fhe 
radiafion  field  quickly  becoming  diffuse,  and  fhaf 
fhere  is  increased  affenuafion  and  scaffering  for 
colder  ice  (Gilberf  and  Schoonmaker  1990,  Voss 
and  Schoonmaker  1992).  Beam  spread  measure- 
menfs,  when  combined  wifh  radiafive  fransfer 
models,  show  promise  as  a  means  of  defermining 
scaffering  coefficienfs  and  phase  funcfions  from 
mulfiply  scaffering  sea  ice. 

MODELS 

If  is  evidenf  from  fhe  observafional  dafa  fhaf 
fhe  optical  properties  of  sea  ice  vary  greafly.  The 
opfical  properfies  vary  spafially  over  scales  of 
only  a  few  mefers  and  fhey  vary  femporally  as 
fhe  ice  cover  metis  in  fhe  summer  and  freezes  in 
fhe  fall.  An  analysis  of  opfical  observations  has 
demonsfrafed  fhaf  fhe  opfical  properfies  of  sea 
ice  are  direcfly  affecfed  by  fhe  sfafe  and  sfrucfure 
of  fhe  ice.  Models  are  essential  in  inferprefing 
observafions  and  in  progressing  from  a  phenom¬ 
enological  collection  of  observafions  fo  a  physi¬ 
cally  based  undersfanding  of  radiafive  fransfer  in 
sea  ice. 

The  variabilify  in  opfical  properfies  also  cre- 
afes  difficulties  in  exfrapolafing  observafions.  In¬ 
dividual  observafions  provide  information  on  fhe 
opfical  properfies  af  a  parficular  locafion  af  a  par¬ 
ticular  time,  buf  for  many  problems  more  general 
information  is  needed  on  how  fhe  opfical  proper¬ 
fies  of  a  region  evolve  wifh  fime.  In  principle,  fhis 
information  can  be  obfained  observafionally,  buf 
for  a  large-scale,  long-ferm  sfudy  fhis  is  nof  prac¬ 
tical.  For  such  sfudies,  models  are  an  essential 
fool. 

Radiafive  fransfer  models  have  been  applied 
fo  a  wide  range  of  problems,  including  esfimaf- 
ing  fhe  absorpfion  of  solar  radiafion  in  sea  ice 
(Maykuf  and  Unfersfeiner  1971),  sfudying  fhe  re- 
lafionship  befween  changes  in  ice  physical  prop¬ 


erfies  and  changes  in  opfical  properfies  (Grenfell 
1983,  1991),  analyzing  fhe  spread  of  a  beam  of 
lighf  as  if  passes  fhrough  ice  (Trodahl  ef  al.  1987), 
investigating  bio-opfical  inferacfions  (Arrigo  ef 
al.  1991),  examining  fhe  fransmission  of  visible 
and  ulfraviolef  lighf  fhrough  sea  ice  (Perovich 
1990,  1991,  1993)  and  assessing  radiafive  inferac¬ 
fions  between  the  atmosphere,  ice  and  ocean  (Jin 
et  al.  1994). 

These  radiative  transfer  models  for  sea  ice 
range  in  complexify  from  a  simple  wavelengfh- 
infegrafed  parameferizafion  of  an  exponenfial 
decay  law  fo  numerically  infricafe  solutions  of 
fhe  radiance  field  in  fhe  ice.  There  are  several 
differenf  models  wifh  a  variefy  of  solufion 
schemes  and  differenf  inpuf  and  oufpuf  param- 
efers;  however,  fhe  same  physics  underlies  all  of 
fhese  models.  They  may  use  differenf  fechniques 
buf  fhey  all  freaf  fhe  basic  physical  properfies  of 
absorpfion  and  scaffering  of  lighf  in  fhe  ice.  Be¬ 
cause  of  fheir  diversify,  fhese  models  all  have 
affribufes  fhaf  endorse  fhem  for  some  applica- 
fions  and  resfricf  fhem  for  ofhers.  A  sampling  of 
sea  ice  radiafive  fransfer  models  is  presenfed  in 
Table  2. 

One  of  fhe  disfinguishing  feafures  of  radiafive 
fransfer  models  is  fhe  number  of  "sfreams"  fhey 
consider.  The  number  of  sfreams  refers  fo  fhe  num¬ 
ber  of  momenfs  from  which  fhe  radiance  is  calcu- 
lafed.  Quife  common  are  fwo-sfream  models, 
where  fhe  upwelling  and  downwelling  irradiances 
are  compufed.  More  sfreams  means  more  angu¬ 
lar  defail  in  fhe  calculafed  radiance  field.  The  cosf 
of  fhis  additional  defail  is  more  complexify  in  fhe 
compufafions  and  offen  a  requiremenf  for  more 
defailed  informafion  on  fhe  opfical  properfies  of 
fhe  ice. 

The  simplesf  sea  ice  radiafive  fransfer  model  is 
fhe  exponenfial  decay  relationship 

F(z,X)  =  (l-aOFoWe"''""  (5) 

where  Fq(X)  is  fhe  incidenf  solar  irradiance.  This 
formulafion  has  fhe  advanfage  of  being  simple 
compufafionally.  However,  fhere  is  an  implicif 
assumpfion  fhaf  fhe  medium  is  infinifely  fhick, 
and  consequenfly,  fhe  exponenfial  decay  law  does 
a  poor  job  of  represenfing  radiafive  fransfer  in 
fhin  ice  (Grenfell  1979). 

Grenfell  (1979)  developed  a  fwo-sfream,  fhree- 
layer  model,  based  on  fhe  work  of  Dunkle  and 
Bevans  (1957),  fhaf  improved  fhe  freafmenf  of 
fhin  ice  wifh  only  a  modes!  increase  in  compufa- 
fional  complexify.  The  general  solufion  for  fhe 
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Table  2.  Summary  of  sea  ice  radiative  transfer  models. 


Model 

Streams 

Spectral 
range  (nm) 

Number 
of  layers 

Solution 

scheme 

Output 

parameters 

Comments 

Grenfell  (1979) 

2 

400-2150 

3 

Analytic 

Fd,  Fu'  ^ 

Examined  thin  ice,  devel¬ 
oped  parameterizations 
for  a  and  F  as  a  function 
of  thickness,  isotropic 
scattering 

Perovich  and  Grenfell  (1982)  14 

400-1000 

2 

DOM*-analytic 

7(6),Fd,F^,a,T 

Anistropic  scattering,  esti¬ 
mate  scattering  param¬ 
eters  from  observations 
of  a  and  T 

Grenfell  (1983) 

16 

350-2750 

1 

DOM-numerical 

l(0),Fd,Fu.a,r 

Detailed  angular  resolution, 
optical  and  physical  pro¬ 
perties  are  related 

Trodahl  et  al.  (1987) 

500,  700 

multiple 

Monte  Carlo 

I  (Q,x) 

Isotropic  and  anisotropic 
scattering,  treats  beam 
spread 

Perovich  (1990, 1993) 

2 

250-1000 

multiple 

Analytic 

Fd/  Fjj,  a,  T 

Ultraviolet  and  visible  wave¬ 
lengths,  computationally 
simple,  easy  ice  charac¬ 
terization,  isotropic  scat¬ 
tering 

Arrigo  et  al.  (1991) 

1 

400-700 

multiple 

Exponential 

Fd,r 

Detailed  treatment  of  im¬ 
pact  of  biogenic  material 
on  light  transmission 

Grenfell  (1992) 

4 

350-2750 

multiple 

DOM-analytic 

Tied  closely  to  ice  physical 
properties,  treats  vertical 
variability  in  ice 

Jin  et  al.  (1994) 

select 

250-4000 

multiple 

DOM-numerical 

1(B),  F^,  7^,  a,  T 

Coupled  atmosphere-ice- 
ocean  radiative  transfer 
absorption  model,  deter¬ 
mines  solar  absorption  in 
each  component 

""Discrete  ordinates  method  (Chandrasekhar  1960) 


up  welling  (F^)  and  downwelling  (Fj)  irradiances 
are 

F(j  (z,  X)  =  A  sinh  (k;^z)  +  B  cosh  (k;^z) 

F^  (z.  A.)  =  C  sinh  (k;^z)  +  B  cosh  (k;^z) 

where  A,  B,  C  and  D  are  determined  from  the 
boundary  conditions.  For  an  optically  thick  me¬ 
dium  (z  ^  oo),  this  solution  converges  to  the  ex¬ 
ponential  decay  law  (eq  5).  The  major  deficiency 
of  the  two-stream  model  is  its  treatment  of  scat¬ 
tering,  in  particular,  the  simplifying  assumption 
of  isotropic  scattering.  An  important  advantage 
of  this  formulation  is  that  it  directly  utilizes  the 
observations  of  light  extinction  in  sea  ice  made 
by  Grenfell  and  Maykut  (1977)  and  Perovich  and 
Grenfell  (1981).  Because  of  this,  only  a  qualitative 
description  of  the  ice  is  needed;  blue  or  white, 
melting  or  cold,  snow-covered  or  bare.  Grenfell 


(1979)  used  this  model  to  investigate  the  depen¬ 
dence  of  albedo,  transmittance,  and  fg  on  thick¬ 
ness  and  ice  t5q)e.  He  then  used  the  results  to 
derive  simple  parameterized  formulae  for  aj  and 
ig  suitable  for  use  in  sea  ice  thermodynamic 
models. 

This  two-stream  formulation  was  expanded 
into  an  n-layer  model  (Perovich  1990)  and  ex¬ 
tended  into  the  ultraviolet  (Perovich  1993).  The 
focus  of  these  studies  was  on  the  spatial  and  tem¬ 
poral  variability  of  reflection,  absorption,  and 
transmission  of  solar  radiation  by  sea  ice.  To  il¬ 
lustrate  the  utility  of  such  models  let  us  examine 
a  particular  problem  of  interest:  the  transmission 
of  visible  and  ultraviolet  light  through  sea  ice  in 
the  Weddell  Sea  during  spring.  Spring  is  the  pe¬ 
riod  when  ozone  depletion  is  the  greatest,  as  is 
the  consequent  increase  in  incident  ultraviolet  ir- 
radiance  and  potential  biological  hazard.  Physi- 
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Figure  17.  Theoretical  estimates  of  ultra¬ 
violet  and  visible  light  transmission 
through  sea  ice  in  the  Weddell  Sea.  Ice 
thickness,  snow  depth  and  physical  prop¬ 
erties  data  are  from  Lange  and  Eicken 
(1991).  Transmitted  UV-B  and  visible 
irradiance  were  computed  using  a  two- 
stream  model  (Perovich  1990, 1993) 


cal  properties  data  from  Lange  and  Eicken  (1991) 
were  used  to  define  fhe  fype  and  fhickness  of  fhe 
ice  and  snow  cover  (Fig.  17a)  along  a  100-m 
fransecf  in  fhe  Weddell  Sea.  Ice  fhicknesses  in  fhis 
area  varied  from  0.2  fo  1.2  m,  while  fhe  snow 
depfh  ranged  from  0.0  fo  0.2  m.  Wifh  fhese  inpuf 
paramefers,  fhe  model  calculafed  esfimafes  of 
fransmiffance  for  fhe  biologically  harmful  UV-B 
irradiance  (280  fo  320  nm)  and  fhe  beneficial  pho- 
fosynfhefically  active  radiafion  (400  fo  700  nm) 
(Fig.  17b).  There  is  fremendous  spafial  variabilify 
in  fhe  UV-B  fransmiffance  over  fhe  100-m  fransecf, 
wifh  values  ranging  over  nearly  fwo  orders  of 
magnifude  from  0.0015  fo  0.09.  The  primary  in¬ 
fluence  on  fransmiffance  is  fhe  snow  depfh,  fol¬ 
lowed  by  fhe  ice  fhickness.  Maximum  fransmif- 
fances  are  associafed  wifh  minimum  snow  depfhs. 
If  is  evidenf  fhaf  fhe  presence  of  an  ice  cover 
causes  a  marked  reduction  in  fransmiffed  lighf 
levels.  This  reducfion  is  greafer  for  fhe  harmful 
UV-B  fhan  for  fhe  beneficial  visible,  implying  fhaf 
sea  ice  may  moderafe  fhe  biological  impacf  of 
enhanced  incidenf  ulfraviolef  irradiance  on  biofa 
living  in  and  under  fhe  ice. 

Models  based  on  fhe  discrefe  ordinafes  mefhod 
(DOM)  of  Chandrasekhar  (1960)  have  been  used 
fo  freaf  scaffering  in  more  defail  and  examine  fhe 
angular  disfribufion  of  radiance.  In  fhe  DOM,  fhe 
phase  funcfion  is  approximafed  by  a  series  of 
Legendre  polynomials  (Liou  1973,  1974,  Mobley 
1994).  The  discrefe  ordinafes  refer  fo  parficular 
angles  af  which  fhe  radiance  is  compufed.  These 
angles  are  nof  arbifrary,  buf  are  defermined  from 


fhe  roofs  of  fhe  Legendre  pol5momial.  In  fhis  for- 
mulafion,  if  is  no  longer  necessary  fo  assume  fhaf 
fhe  radiance  field  is  diffuse  and  fhaf  fhe  phase 
funcfion  is  isofropic.  However,  fhese  models,  par¬ 
ticularly  for  larger  numbers  of  sfreams,  are  sig- 
nificanfly  more  complex  compufafionally. 

Perovich  and  Grenfell  (1982)  developed  a  fwo- 
layer,  four-sfream  model  (radiances  af  fwo  up¬ 
ward  and  fwo  downward  angles)  and  applied  if 
fo  invesfigafe  fhe  effecfs  of  ice  fhickness,  and  fhe 
influence  of  direcf  vs.  diffuse  incidenf  solar  ra¬ 
diafion,  on  specfral  albedo  and  fransmiffance. 
Using  experimenfally  defermined  phase  functions 
fhey  found  fhaf  single  scaffering  albedos  (Oo)  for 
young  ice  were  high:  from  0.95  for  warm  melting 
young  ice  fo  0.9997  for  young  ice  below  fhe  eu- 
fecfic  poinf. 

Grenfell  (1983,  1991)  developed  a  single-layer, 
16-sfream  model  and  a  multilayer,  four-sfream 
model  fo  explore  relationships  befween  ice  physi¬ 
cal  properties  and  ice  optical  properties.  The  four- 
sfream  (Grenfell  1991)  model  significanfly  ex- 
fended  fhe  work  of  Perovich  and  Grenfell  (1982) 
by  including  vertically  varying  ice  properties.  The 
single-layer,  16-sfream  model  (Grenfell  1983)  gen- 
erafed  a  more  defailed  angular  description  of  ra¬ 
diance,  beffer  represenfed  fhe  phase  funcfion,  and 
improved  fhe  freafmenf  of  refracfion  af  fhe  air-ice 
inferface  for  a  homogeneous  ice  cover.  This  model 
was  used  fo  direcfly  link  fhe  physical  properties 
of  fhe  ice,  such  as  fhe  inclusion  size  disfribufions 
of  air  bubbles  and  brine  pockefs,  fo  radiative  frans- 
fer  in  fhe  ice.  The  absorpfion  and  scaffering  coef- 
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Figure  18.  Calculated  estimates  of  spectral  albedo  as  a  function  of 
ice  density  and  growth  rate  (from  Grenfell  1983).  The  ice  was  3 
m  thick.  The  air  volume  was  zero  for  the  growth  rate  simulation. 
Figure  18a  shows  albedo  as  a  function  of  ice  density  (p):  a)  p  = 
0.86  g  cm~^,  b)  p=0.88  g  cm~^,  c)  p  =  0.90  g  cm~^,  d)  p  =  0.91  g 
cm~^  and  e)  p=  0.94  g  cm~^.  Figure  18  shows  albedo  as  a  func¬ 
tion  of  growth  rate  (f)  for  a)  f  =  8  x  10~^  cms~^,b)f=  4  x  10~^ 
cm  s~^,  c)f=2x  ICr^  cm  s~^  and  d)  f  =  8  xlCr^  cm  s~^ . 


ficients  depended  explicitly  on  the  amount  and 
size  distribution  of  air  bubbles  and  brine  pockets. 
These  values  in  turn  depended  on  the  ice  growth 
conditions,  thermal  history,  temperature,  salinity 
and  density  With  this  formulation,  it  was  pos¬ 
sible  to  theoretically  explore  the  impact  of  growfh 
conditions  and  fhermal  hisfory  on  specfral  albe¬ 
dos  and  extinction  coefficienfs.  For  example.  Fig¬ 
ure  18  shows  calculafed  esfimafes  of  specfral  al¬ 
bedo  for  differenf  ice  densifies  and  ice  growfh 
rafes.  The  ice  was  3  m  fhick  in  fhese  cases.  The 
large  impacf  of  air  bubbles  on  scattering  and  ice 
opfical  properfies  is  demonsfrafed  in  Figure  18a. 
There  is  an  increase  in  albedo  as  fhe  ice  densify 
decreases  and  fhe  number  of  air  bubbles  increases. 
This  increase  is  mosf  pronounced  af  470  nm,  where 
absorption  is  smallesf.  The  albedo  af  470  nm  was 
abouf  0.57  for  bubble-free  ice  (p  =  0.94)  and  in¬ 
creased  fo  0.84  for  bubbly  ice  wifh  an  air  volume 
of  8%  (p  =  0.86).  Calculations  also  indicafed  fhaf 
fasfer  growfh  rafes  resulf  in  larger  albedos  (Fig. 
18b).  For  fhese  calculafions  fhe  air  volume  was 
assumed  fo  be  zero,  so  changes  in  albedo  resulfed 
from  changes  in  fhe  plafelef  spacing  and  fhe  num¬ 
ber  of  brine  inclusions.  Fasfer  grown  ice  has 


smaller  plafelefs,  higher  salinify,  and  more 
brine  inclusions  (Weeks  and  Ackley  1982). 
This  is  consisfenf  wifh  our  premise  fhaf 
more  inclusions  means  more  scattering 
and  higher  albedos. 

The  Monfe  Carlo  mefhod  is  anofher 
approach  fo  radiafive  fransfer  modeling. 
As  fhe  name  implies,  Monfe  Carlo  models 
fake  a  sfafisfical  approach  fo  solving  fhe 
equation  of  radiafive  fransfer  (eq  2).  In 
shorf,  fhe  absorpfion  coefficienf,  fhe  scaf- 
fering  coefficienf  and  fhe  phase  funcfion 
are  fransformed  info  fhe  probabilify  fhaf 
over  a  given  disfance  a  photon  is  absorbed 
or  scaffered,  and  if  scaffered,  in  whaf  di- 
recfion.  Wifh  fhese  probabilifies  known, 
enormous  numbers  of  photons  are  numeri¬ 
cally  "shof"  info  fhe  medium.  The  fafe  of 
each  phofon  is  decided  by  fhe  roll  of  fhe 
dice,  or  more  precisely,  fhe  whim  of  fhe 
random  number  generafor.  Radiafive 
fransfer  in  fhe  medium  is  described  by 
fhe  cumulafive  resulf  of  all  fhe  phofons. 
Because  of  fhe  large  number  of  phofons 
needed,  Monfe  Carlo  models  are  very  in- 
efficienf  compufafionally.  They  are,  how¬ 
ever,  simple  concepfually,  simple  fo  pro¬ 
gram,  and  widely  applicable  (Mobley 
1994).  This  mefhod  is  particularly  well 
suited  for  complex  geomefries  or  boundary  con¬ 
ditions,  where  ofher  solufions  fo  fhe  equafion  of 
radiafive  fransfer  are  difficulf  or  impossible. 
Trodahl  ef  al.  (1987)  and  Trodahl  and  Buckley 
(1989)  effectively  used  Monfe  Carlo  solufions  in 
beam  spread  sfudies,  bofh  fo  model  observations 
and  infer  informafion  on  fhe  scattering  proper¬ 
fies  of  sea  ice.  They  found  fhaf  scaffering  in  fhe 
surface  layer  of  fhe  ice  was  greafer  fhan  in  fhe 
inferior  and  fhaf  fhe  scaffering  was  anisofropic. 

An  excifing  new  modeling  developmenf  has 
been  fhe  inclusion  of  biological  effecfs  in  sea  ice 
opfical  models.  Sea  ice  is  fhe  habifaf  of  a  rich 
microbial  communify  (Palmisano  and  Sullivan 
1983,  Garrison  ef  al.  1986).  Ice  biofa  bofh  affecf 
and  are  affecfed  by  fhe  specfral  irradiance  wifhin 
fhe  ice.  Arrigo  ef  al.  (1991)  developed  a  bio-opfi- 
cal  model  fo  invesfigafe  fhe  inferdependence  be- 
fween  biology  and  fransmiffed  lighf.  They  used  a 
simple  exponential  decay  law  fo  model  irradi¬ 
ance  wifhin  fhe  ice,  buf  fhey  coupled  fhis  wifh  a 
sophisficafed  freafmenf  of  fhe  exfincfion  coeffi¬ 
cienfs  (k).  They  formulafed  polynomial  relation¬ 
ships  defining  specfral  exfincfion  coefficienfs  for 
dry  snow,  wef  snow,  congelafion  ice,  plafelef  ice. 


18 


Figure  19.  Seasonal  changes  in  underice  spectral  irra- 
diance  calculated  using  a  bio-optical  model  (Arrigo  et 
al.  1991).  Curves  are  predicted  spectral  transmitted 
irradiance  at  noon  on  a)  7  October  1984,  b)  13  Novem¬ 
ber  1984  and  c)  5  December  1984. 

ice  cooler  than  the  eutectic  point,  and  as  a  func¬ 
tion  of  brine  volume.  Mosf  imporfanfly,  fhey  also 
derived  relafionships  for  fhe  exfincfion  confribu- 
fions  from  absorpfion  due  fo  microalgae  and  de- 
frifus.  Wifh  fhis  model  if  is  possible  fo  examine 
fhe  impacf  of  biogenic  maferial  on  fransmiffed 
specfral  irradiance  and  fo  invesfigafe  femporal 
changes. 

Combining  field  observafions  wifh  model  cal- 
culafions,  Arrigo  ef  al.  (1991)  were  able  fo  com- 
pufe  a  fime  series  of  fransmiffed  specfral  irradi- 
ances  (Fig.  19).  The  calculafions  were  done  for 
bare  ice  in  McMurdo  Sound,  Anfarcfica,  roughly 
1.7-1. 8  m  fhick,  between  7  October  and  5  Decem¬ 
ber.  During  fhis  period  fhere  was  a  consfanf  in¬ 
crease  in  fhe  amounf  of  microalgae  and  defrifus. 
On  7  Ocfober,  levels  of  microalgae  were  low  and 
fhere  was  no  defrifus  presenf,  so  lighf  losses  were 
primarily  due  fo  exfincfion  by  fhe  sea  ice.  By  13 
November  fhe  spring  bloom  had  produced  sig- 
nificanf  amounf s  of  algae  and  defrifus.  The  pres¬ 
ence  of  fhis  biogenic  maferial  resulfed  in  an  over¬ 
all  reducfion,  and  a  change  in  fhe  specfral  shape, 
of  fhe  fransmiffed  irradiance.  The  disfincf  spec¬ 


fral  shape  of  fhe  fransmiffed  irradiance  is  charac- 
ferisfic  of  ice  wifh  biogenic  maferial.  Algae  and 
defrifus  levels  confinued  fo  increase  fhrough  5 
December,  causing  a  furfher  reducfion  in  frans- 
miffance. 

A  model  was  recenfly  developed  fo  examine 
radiafive  fransfer  in  a  coupled  afmosphere-ice- 
ocean  system  (Jin  ef  al.  1994).  The  model  is  a 
multilayer  and  mulfisfream  formulation  based  on 
fhe  discrete  ordinafes  mefhod.  Radiafive  fransfer 
wifhin  fhe  entire  afmosphere-ice-ocean  system  is 
defermined  based  on  a  description  of  physical 
properties  of  fhe  afmosphere,  ice  and  ocean  from 
which  fhe  opfical  properfies  are  derived.  The 
model  compufes  fhe  disfribufion  and  absorpfion 
of  solar  radiation  in  fhe  afmosphere,  ice  and  ocean. 
Resulfs  indicate  fhaf  sea  ice  has  a  sfrong  influ¬ 
ence  on  fhe  disfribufion  of  solar  radiation  in  fhe 
sysfem  (Jin  ef  al.  1994).  Such  models  provide  a 
promising  fool  for  invesfigafing  afmosphere-ice- 
ocean  radiafive  feedbacks. 


SUMMARY  AND  CURRENT 
AREAS  OF  INTEREST 

By  now  fhe  reader  is  no  doubf  aware  fhaf  fhe 
opfical  properfies  of  sea  ice  are  variable  and  com¬ 
plex.  The  reader  is  also  aware  fhaf  much  of  fhis 
complexify  is  comprehensible.  While  many  of  fhe 
defails  still  need  fo  be  defermined,  we  do  have  a 
qualifafive  undersfanding  of  sea  ice  opfical  prop¬ 
erfies  and  fheir  variabilify.  This  undersfanding  is 
based  on  a  few  fundamenfal  principles.  Changes 
in  such  opfical  properfies  as  fhe  albedo,  reflec- 
fance,  fransmiffance,  and  exfincfion  coefficienf  are 
direcfly  fied  fo  changes  in  fhe  sfafe  and  sfrucfure 
of  fhe  ice.  Physical  changes  in  fhe  ice  which  en¬ 
hance  scaffering,  such  as  fhe  formafion  of  air 
bubbles  due  fo  brine  drainage,  resulf  in  larger 
albedos  and  exfincfion  coefficienfs.  Radiafive 
fransfer  in  sea  ice  is  a  combinafion  of  absorpfion 
and  scaffering.  Differences  in  fhe  magnifude  of 
fhese  opfical  properfies  are  due  primarily  fo  dif¬ 
ferences  in  scaffering.  Specfral  variations  are 
mainly  a  resulf  of  absorpfion. 

In  addition  fo  fhese  general  principles,  fhere 
are  also  several  specific  commenfs  fhaf  can  be 
made  regarding  fhe  sea  ice  opfical  properfies. 
Specfral  absorpfion  coefficienfs  for  ice  are  well 
known,  however,  represenfafive  values  for  brine 
are  less  cerfain.  Absorpfion  by  algae  and  parficu- 
lafes  is  also  imporfanf  and  needs  furfher  invesfi- 
gafion.  More  work,  bofh  experimenfal  and  fheo- 
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retical,  needs  to  be  done  investigating  scattering 
in  sea  ice.  The  albedo  is  quite  sensitive  to  the 
surface  state.  If  fhe  ice  has  an  appreciable  snow 
cover,  visible  wavelengfh  albedos  are  above  90% 
and  little  lighf  is  fransmiffed  fo  fhe  ocean.  In  very 
cold  ice  (T<-24°C),  hydrohalife  precipifafes,  caus¬ 
ing  a  sharp  increase  in  albedo  and  exfincfion  co- 
efficienf  fo  values  comparable  fo  snow.  There  is  a 
less  pronounced,  buf  sfill  pofenfially  significanf, 
effecf  af  femperafures  below  -8°C  where  mirabilife 
precipifafes.  The  presence  of  liquid  wafer  on  fhe 
surface  causes  a  decrease  in  albedo,  which  is  more 
pronounced  af  longer  wavelengfhs.  If  fhe  surface 
drains,  fhe  brine  pockefs  become  air  bubbles,  re- 
sulfing  in  more  scattering  and  an  increase  in  al¬ 
bedo  and  exfincfion  coefficienf.  Ice  fhaf  is  grown 
fasfer  has  more  plafelefs  and  more  brine  inclu¬ 
sions,  and  consequenfly  large  albedos  and  ex¬ 
fincfion  coefficienfs.  The  opfical  properfies  de¬ 
pend  nof  only  on  fhe  volume  of  brine  or  air,  buf 
on  how  fhaf  brine  or  air  is  disfribufed. 

Sea  ice  opfical  properfies  is  currenfly  a  research 
area  of  considerable  inferesf  and  acfivify.  Even 
fhough  much  has  been  learned  abouf  fhe  opfical 
properfies  of  sea  ice,  fhere  are  sfill  numerous  im- 
porfanf  and  infriguing  problems  exfanf.  A  major 
goal  is  quanfifying  relafionships  befween  fhe 
physical  and  fhe  opfical  properfies.  Achieving  fhis 
goal  enfails  nof  only  a  better  undersfanding  of 
fhe  opfical  properfies,  buf  a  beffer  undersfanding 
of  fhe  physical  properfies.  Because  of  fhe  pofen- 
fial  climafological  impacf  of  ice-albedo  feedback, 
one  area  of  parficular  concern  is  defermining  how 
fhe  changes  in  fhe  physical  sfafe  of  fhe  ice  during 
fhe  summer  melf  season  affecf  fhe  albedo  of  fhe 
ice  cover. 

An  improved  undersfanding  of  scattering  in 
sea  ice  is  needed.  This  can  be  addressed  fhrough 
laborafory  sfudies  of  fhe  scattering  properfies  of 
small  sea  ice  samples  (Miller  ef  al.  1994)  and 
fhrough  field  sfudies  invesfigafing  fhe  spread  of 
a  collimafed  beam  of  lighf  in  ice  (Longacre  and 
Landry  1994).  Anofher  approach  fo  esfimafing  fhe 
scaffering  properfies  of  sea  ice  is  fo  use  Mie  fheory 
(Bohren  and  Huffman  1983).  Asfafisfical  descrip¬ 
tion  of  fhe  ice  microsfrucfure  is  needed  for  fhis 
approach,  including  defailed  informafion  on  fhe 
inclusion  size  disfribufions  for  fhe  air  bubbles 
and  brine  pockefs  (Perovich  and  Gow  1991).  Little 
is  known  regarding  fhese  size  disfribufions  and 
how  fhey  vary  wifh  ice  physical  properfies  such 
as  brine  volume,  densify  and  growfh  rafe. 

In  fhe  pasf  fhere  has  been  an  abundance  of 
albedo  measuremenfs,  buf  few  observafions  of 


fransmiffed  lighf.  This  deficiency  has  impeded 
radiafive  fransfer  modeling  efforfs,  ice  heaf  bal¬ 
ance  sfudies,  and  bio-opfical  invesfigafions.  This 
is  begirming  fo  change  as  advancing  fechnology 
leads  fo  improved  insfrumenfafion  and  irmova- 
five  new  approaches  fo  measuring  lighf  in  and 
under  fhe  ice  are  developed.  New  sensors  make 
if  possible  fo  measure  defailed  specfral  fransmif- 
fances  even  under  fhick  snow-covered  ice.  Fiber 
opfic  probes  can  be  frozen  in  fhe  ice  fo  measure 
fhe  radiance  disfribufion  wifhin  fhe  ice.  Powerful 
fechniques  are  being  applied  fo  measure  in-sifu 
profiles  of  fransmiffed  irradiance,  beam  spread, 
and  diffuse  affenuafion  coefficienf.”^ 

Many  pressing  issues  concerning  sea  ice  opfi¬ 
cal  properfies  can  only  be  addressed  fhrough  in- 
ferdisciplinary  sfudies.  A  combined  efforf  is 
needed  fo  examine  such  issues  as  assessing  ice- 
albedo  feedback,  ascerfaining  fhe  impacf  of  en¬ 
hanced  incidenf  levels  of  ulfraviolef  irradiance 
on  biofa  living  in  or  under  fhe  ice,  and  using 
safellife-measured  microwave  signafures  as  a 
proxy  for  large-scale  ice  albedo.  Recenf  experi- 
menfal  programs  have  recognized  fhis  and  have 
emphasized  acquiring  a  comprehensive  dafa  sef, 
including  informafion  on  fhe  ice  sfafe  and  sfruc- 
fure,  biofa,  parficulafes  and  microwave  signafures, 
as  well  as  complefe  opfical  measuremenfs. 

Anofher  approach  fo  fhese  problems  is  fhrough 
modeling,  in  parficular  fhrough  fhe  infegrafion 
of  models.  As  fhe  previous  secfion  demonsfrafed, 
fhere  are  several  good  radiafive  fransfer  models 
for  sea  ice  fhaf  include  informafion  on  fhe  physi¬ 
cal  properfies  of  fhe  ice  (Grenfell  1983, 1991,  Jin  ef 
al.  1994).  There  are  also  models  fhaf  freaf  fhe 
physical  properfies  of  sea  ice  during  fhe  firsf  year 
of  growfh  (Gox  and  Weeks  1988,  Wade  and  Weeks 
in  press).  Thermodynamic  sea  ice  models  include, 
fypically  in  a  parameferized  fashion,  fhe  reflec- 
fion,  absorption  and  fransmission  of  solar  radia¬ 
tion.  The  effecfs  of  biogenic  maferial  on  fransmif¬ 
fed  irradiance  can  be  considered  (Arrigo  ef  al. 
1991),  and  fhere  has  been  progress  fowards  de¬ 
veloping  a  frue  bio-opfical  model  where  fhe  infri- 
cafe  inferplay  befween  fhe  lighf  levels  in  and  un¬ 
der  fhe  ice  and  fhe  amounf  of  biological  acfivify 
can  be  fully  explored  (Arrigo  ef  al.  1993). 

General,  comprehensive,  inferdisciplinary 
models  are  needed  models  fhaf  couple  fhe  ice 


*  Personal  communication  with  S.  Pegau,  College  of 
Oceanic  and  Atmospheric  Sciences,  Oregon  State  Uni¬ 
versity. 
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physical  properties,  optical  properties,  biological 
properties  and  thermodynamics,  and  that  link  the 
ice  to  the  atmosphere  and  ocean.  Models  where 
changes  in  ice  temperature  cause  changes  in  the 
physical  properties  of  the  ice,  which  in  turn  im¬ 
pact  the  ice  optical  properties  and  thereby  the 
physical  properties  and  the  biological  activity 
within  the  ice.  Models  where  changes  in  the  ice 
are  coupled  to  energy  exchange  with  the  atmo¬ 
sphere  and  ocean.  Developing  such  models  is  a 
daunting  task,  but  a  task  with  substantial  rewards. 
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APPENDIX  A:  LIST  OF  SYMBOLS 
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radiance  of  the  direct  beam  component  of  the  incident  radiation  field 

growth  rate 

irradiance 

downwelling  irradiance 
upwelling  irradiance 
radiance 

fraction  of  incident  irradiance  transmitted  through  the  top  0.1  m  of  the  ice 

reflected  radiance 

absorption  coefficient  of  brine 

absorption  coefficient  of  ice 

absorption  coefficient  of  sea  ice 

real  index  of  refraction 

phase  function 

bidirectional  reflectance  distribution  function  (BRDF) 

normalized  reflected  radiance  at  nadir 

source  function 

transmittance 

horizontal  position 

depth  within  the  medium 

albedo 

wavelength-integrated,  or  total,  albedo 
azimuth  angle 
solar  azimuth  angle 
extinction  coefficient 

wavelength-integrated,  or  total,  extinction  coefficients 

wavelength 

volume  fraction  of  ice 

volume  fraction  of  brine 

zenith  angle  (0  pointing  downward,  n  pointing  upward) 

solar  zenith  angle 

cosine  of  the  zenith  angle,  6 

density 

scattering  coefficient 
nondimensional  optical  depth 
single  scattering  albedo 
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